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Abstract 

Supersymmetry is an attractive extension of the standard model of particle 
physics. It associates to every bosonic degree of freedom a fermionic one and 
vice versa. Supersymmetry unifies the coupling constants of the electromag- 
netic, weak and strong forces at a high scale and provides a candidate for 
the elusive dark matter. Supersymmetry could be discovered at the LHC, 
the proton-proton collider at CERN which has started operations in 2008. 
The LHC is foreseen to have a center-of-mass energy of 14 TeV, opening 
up a new mass range to be explored to search for supersymmetric particles 
with the ATLAS and CMS experiments. The development and production 
of electronics for these detectors has been a challenge, e.g. for the readout 
board for the electromagnetic calorimeter. Reconstructing the physics events 
with the best precision, in particular the reconstruction and identification of 
electrons and photons in the large QCD background has been prepared in 
extensive test beam studies and Monte Carlo simulations. If the Higgs boson 
and supersymmetry are discovered, the properties of the (s)particles such as 
the masses, branching ratios must be measured precisely, either at the LHC 
or at a future e^e~ linear collider. The SFitter project aims to determine 
the underlying theoretical model parameters from the correlated experimen- 
tal measurements including theoretical errors. The methods are applied to 
the extraction of the fundamental parameters of supersymmetry as well the 
measurement of the Higgs boson couplings at the LHC. The extrapolation of 
the supersymmetric parameters from the weak scale to the Grand Unification 
Scale could provide the basis towards the inclusion of gravity. 



Therefore, since brevity is the soul of wit, 

And tediousness the limbs and outward flourishes, 

I will be brief. 



Hamlet, Act 2 Scene 2, William Shakespeare 
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Chapter 1 
Introduction 



The Standard Model of elementary particle physics [H El [3] describes matter 
and its interactions with an unprecedented degree of precision. Matter is 
built of fermions (quarks and leptons) and their interactions are mediated via 
bosons. The photon (7) is responsible for the electromagnetic interactions, 
the charged vector bosons (W^) and the Z boson are the carriers of the weak 
force. The gluons mediate the strong interaction. Masses are generated via 
the Higgs mechanism, leading to an additional neutral scalar particle, the 
Higgs boson (H) m El El [a E]. 

Supersymmetry is the next logical step after gauge theories [9]. In su- 
persymmetric theories a fermionic degree of freedom is associated to every 
bosonic degree of freedom and vice-versa [101 El 1121 lEl [H] (and references 
therein). Ultra-violet completions of the Standard Model lead to quantum 
corrections on the Higgs boson mass of the order of the new physics scale, 
larger than the Higgs boson itself, whereas in supersymmetry these large cor- 
rections cancel. Supersymmetry can solve this hierarchy problem in a natural 
way. In the Standard Model the coupling constants of the strong, weak and 
electromagnetic interactions, extrapolated via renormalization group equa- 
tions (RGEs), do not meet in one point. Supersymmetry modifies this behav- 
ior to unify the constants at about 10^^ GeV (grand unification scale, GUT 
scale) [ini [m [17] . Additionally supersymmetry is attractive as it provides a 
candidate for the elusive dark matter. 

To construct supersymmetric models such as the minimal supersym- 
metric extension of the Standard Model (MSSM) or minimal Supergravity 
(mSUGRA), the Higgs sector has to be extended to give masses to up and 
down type quarks The two Higgs doublets are necessary: to avoid anomalies 



and large flavor changing neutral currents (FCNC). In the Standard Model 
the Higgs field and its complex conjugate field are used to give masses to up 
and down type quarks, while in a supersymmetric theory the complex conju- 
gate field cannot be used, necessitating the second doublet. The Higgs sector 
of the MSSM contains five Higgs bosons: h, H, A and H^. The lightest Higgs 
boson must be light with a mass of less than 140 GeV/c^. Additionally the 
sfermions, the partners of the fermions, are predicted. The supersymmetric 
partners of the neutral gauge and Higgs bosons mix to form the neutralinos. 
The charginos are mixtures of the supersymmetric partners of the charged 
Higgs and gauge bosons. 
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Figure 1.1: The expected and ob- 
served lower limit on the Standard 
Model Higgs boson mass obtained by 
the LEP collaborations is shown [18] . 





Tevatron Run II Preliminary, L=0.g-4.2 fb'^ 




LEP Exclusion Tevatron 




f.. Pv^hicinri 


10 


Expected 

Observed ^ 

■■ >1n Expected 

2" Ex poet fd 


1 


.....---•'^^-^ .......:::?^ 


SM 

7 "rn ?=«?,,,,- 



100 110 120 130 140 150 160 170 180 190 200 

mn(GeV/c^) 

Figure 1.2: The expected and ob- 
served lower limit on the Standard 
Model Higgs boson mass obtained by 
the CDF and Do collaborations is 
shown 1191. 



The most precise measurements of parameters of the Standard Model are 
obtained from e^e~ colliders such as LEP and SLC, as well as the proton- 
antiproton (pp) collider TeVatron at Fermilab. ALEPH, DELPHI, L3 and 
OPAL took data from 1989 to 2000 at LEP, operating at center-of-mass en- 
ergies from 91 GeV to 209 GeV. The SLD detector operated at SLAC's e+e" 
SLC at a center-of-mass energy of about 91 GeV. CDF and DO are taking 
data at the TeVatron, operating at a center-of-mass energy of 1.96 TeV. 
In 2009 the proton-proton (pp) collider LHC at CERN will effectively start 



operations. The center-of-inass energy is foreseen to be 7 TeV which will be 
increased to 14 TeV in the coming years. ATLAS and CMS are the two multi- 
purpose detectors located at the LHC They are optimized for the search of 
the Higgs boson(s) and new physics. The e"'"e~ Linear Collider project (ILC) 
is being prepared for a center-of-mass energy from 500 GeV to 1 TeV. A 
decision on whether the ILC will be built is expected after the analysis of the 
first LHC data. R&D is being performed for a linear collider (CLIC) which 
may be capable of running at 3 TeV. 

The search for the Standard Model Higgs boson at LEP has led to a 
limit of 114.4 GeV/c^ [18] as shown in Figure II. ![ slightly smaller than the 
expected limit of 115.3 GeV/c^. Time will tell whether the events observed 
at LEP around this mass will turn out to have been the first indication of 
the Higgs boson or not. 

The TeVatron is particularly sensitive to the Standard Model Higgs boson 
for masses of the order 150-180 GeV/c^. In this mass region the decay of the 
Higgs boson to a pair of W^ bosons is dominant. The combination of several 
channels from both CDF and DO has lead to the first new limit beyond the 
LEP bound: a Standard Model Higgs boson with a mass of 160 to 170 GeV/c^ 
is excluded at 95% confidence level [HI |20] as shown in Figure 11.21 

The precise measurement of the electroweak parameters at LEP and SLC 
are discussed in detail in [21]. The most precise individual measurements 
of the W^ boson mass are from CDF [2^ and DO [23] which have reached 
a precision of 48 MeV/c^ and 43 MeV/c^ respectively. The combination of 
all W^ mass measurements leads to 80.399 ± 0.023 GeV/c^. The top quark 
mass has now been measured with a precision better than percent as shown 
in [21] . The predictions of the electroweak observables depend via quantum 
corrections on the mass of the Higgs boson, especially the W^ boson and 
top quark masses. Their sensitivity is illustrated in Figure 11.41 where the 
solid line denotes the experimental measurements of the W^ boson and top 
quark masses and the straight lines different Higgs boson mass hypotheses. 
Combining all measurements in a fit the Higgs boson mass can be determined 
to be 84l26GeV/c^. This leads to an upper limit of 157 GeV/c^ on the Higgs 
boson mass at 95% confidence level. Thus the data prefer a light Higgs boson 
mass. 

In addition to the Standard Model production mode, albeit with modified 
couplings, the associated production of Higgs bosons, such as the h produced 
in association with the A, can also be used in the search for the MSSM Higgs 
bosons. Limits on neutral MSSM Higgs bosons have been determined to be 
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Figure 1.3: The result of the elec- 
troweak fit is shown as function of 
the Higgs boson mass. The yel- 
low shaded area shows the region ex- 
cluded by the direct search. 
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Figure 1.4: The measurement of the 
W^ boson mass as function of the 
top quark mass is shown (solid line). 
The dotted line shows the indirect 
constraints on the W"*^ boson mass 
and top quark mass. The straight 
lines show different Higgs mass hy- 
potheses (update of [2T]). 



be 92.8 GeV/c^ for the h and 93.4 GeV/c^ for the A boson [25]. The searches 
for the neutral MSSM Higgs bosons can also be interpreted as limit on the 
ratio of the vacuum expectation values (vevs) of the two Higgs doublets tan (3. 
tan/3 between 0.7 and 2 are excluded as shown in Figure [T75l Some variation 
of the excluded region is possible as it depends on the top quark mass and 
the parameters of the supersymmetric model. 

As the cross section for the pair production of scalars is proportional to 
(3^, where (3 = y(l — 4m^±/s), the mass limit on the charged Higgs bosons 
from the combination of the LEP experiments has not reached the kinematic 
limit. Charged Higgs boson masses of less than 78.6 GeV/c^ are excluded at 
95% confidence level [26] . 

LEP has set mass limits on selectrons and muons (cross section ~ /3^) at 
95 GeV/c^ and slightly weaker limits on the stau slepton of 87 GeV/c^ [27] . 
Charginos with masses of up to 103 GeV/c^ (cross section ~ (3), i.e. es- 
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Figure 1.5: The region excluded by 
LEP and by theory in the plane tan /3 
versus mass of the lightest Higgs bo- 
son 
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Figure 1.6: The regions excluded by 
the TeVatron and LEP in the plane 
of the squark mass versus the gluino 
mass are shown [301. 



sentially up to the kinematic limit, have been excluded [28]. Combining all 
searches a lower limit on the mass of the lightest neutralino of 50 GeV/c^ 
has been obtained in mSUGRA [29] . 

The TeVatron is on the forefront of the search for supersymmetric colored 
particles, the squarks and gluinos. In the absence of a signal, lower limits 
on the mass of the sparticles have been set. These depend on the specific 
scenario chosen. As an example, for a given typical choice of parameters, 
squarks of mass less than 392 GeV/c^ have been excluded and a limit of 
327 GeV/c^ on the mass of the gluino has been obtained. For equal squark 
and gluino masses a lower limit of 408 GeV/c^ has been determined [30] (and 
references therein) as shown in Figure 11.61 

Several aspects of the search for supersymmetry and its measurement will 
be discussed in the following. In chapter |2l the phenomenology of the most 
important supersymmetric models will be reviewed. Chapter |3] contains a 
brief discussion of experimental aspects of the search for supersymmetry. In 
particular a part of the Front End Electronics of the ATLAS electromagnetic 
calorimeter as well as the performance of the reconstruction of electrons and 
photons are described in more detail, based on Refs. [311 [321 [33]. In chap- 
ter [H the determination of supersymmetric parameters from experimental 



measurements as well as the determination of the Higgs boson couplings at 
the LHC will be discussed. The results summarized in this chapter are ex- 
plained in more detail in Refs. [Ml [351 EH] . In the last chapter the results 
will be summarized and an outlook for the future will be given. 



Chapter 2 

Phenomenology of 
Supersymmetry 



To construct a supersymmetric extension of the Standard Model, such as the 
MSSM and mSUGRA, an extended Higgs sector is necessary in order to give 
masses to up and down type quarks as well as the addition of partners of 
the standard model particle. Further extensions, such as the MSSM with an 
additional singlet (NMSSM), have an even richer spectrum of Higgs bosons 
and possibly new signatures with respect to the MSSM [37]. In this chapter 
the discussion is restricted to the MSSM and mSUGRA. The parameters of 
the MSSM are defined at the electroweak scale, whereas the parameters of 
mSUGRA are defined at the GUT scale. R-parity will be conserved in the 
following, therefore supersymmetric particles are produced in pairs, (cascade- 
) decaying to the lightest supersymmetric particle (LSP). The LSP is stable, 
neutral and weakly interacting, leading to the characteristic signature for 
supersymmetry at a collider, missing energy. 

The slope of the gauge couplings' running is modified at a scale of about 
1 TeV by Supersymmetry, leading to a unification of the couplings at about 
10^^ GeV. As this scale is close to the Planck scale (10^® GeV), Supersym- 
metry could be a telescope to physics at this scale. A side effect of this 
unification at a higher scale compared to the Standard Model is that the 
proton can have a longer lifetime. 

The Standard Model accounts for only 4% of the energy of the universe. 
24% are due to dark matter. The LSP of supersymmetry is prime candidate 
for dark matter as it is neutral and weakly interacting. 

From an experimental point of view it is intriguing to note that the Stan- 
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dard Model electroweak fits favor a light Higgs boson with a mass less than 
157 GeV/c^. Supersymmetry, as realized in the MSSM, predicts a lightest 
Higgs boson which should have a mass less than about 140 GeV/c^. 

Ultra-violet completions of the Standard Model lead to quantum correc- 
tions on the mass of the Higgs boson of the order of this new physics scale. 
This unnatural behavior is corrected in supersymmetry where an equal num- 
ber of fermions and bosons leads to a cancellation of the corrections. In 
order to preserve the natural cancellation in the Higgs sector the breaking 
of Supersymmetry, i.e. hfting the mass degeneracy of standard and super- 
symmetric particles, cannot be arbitrarily large, but must be at most of the 
order of TeV. 

These arguments motivate the search for Supersymmetry. In the following 
sections a few basic considerations of supersymmetry and the phenomenology 
of two models, the MSSM and mSUGRA, will be discussed. 

2.1 Supersymmetry 

Supersymmetry was presented first in [TOl [TTl [12] . The construction of viable 
models was pioneered by Fayet and Iliopoulos [13] 

At the Planck scale Mp = 2.4 ■ 10^^ GeV/c^, gravitational effects are im- 
portant. The corrections to the Higgs boson mass in the Standard Model 
would have to be fine-tuned to one part in 10^^ to give a Higgs mass of the 
order of the electroweak scale if the Standard Model were to be valid up to 
this scale. As supersymmetric models contain two scalars for every fermion, 
the quadratic divergences cancel and only the logarithmic parts remain. The 
fermion masses, including radiative corrections, are logarithmically diver- 
gent. Including additionally soft supersymmetry breaking, qualitatively the 
corrections to the Higgs boson mass then take on the following form (nigoft 
is the mass splitting between the fermion and the scalars) 



Am^ = m^^ft 



-^ ln(A/msoft) + 



1677 



(2.1) 



If A ~ 1 and A the Planck scale, for a soft supersymmetry breaking mass of 
1 TeV/c^, the correction to the Higgs boson mass will be about 500 GeV/c^. 
Since the correction increases about linearly with nisoft, therefore the masses 
of at least some of the supersymmetric particles must be less than about 
1 TeV/c2. 

10 



Supersymmetry is a spacetime symmetry. The Haag-Lopuzanski-Sohnius 
extension of the Coleman-Mandula [52] theorem restricts the forms of such 
symmetries in interacting quantum field theories. In practice for extensions 
of the Standard Model, the generators Q, Q^ have to satisfy an algebra of 
the form [14] : 

{Q,Q} = {Q^QU = (2.2) 

where Q generates the transformation of a fermion to a boson and vice versa, 



and P^ is the momentum generator of spacetime translations. As in [H] the 
spinor indices have been omitted. 

Supermultiplets are introduced in which a boson and a fermion must be 
present. Scalars and fermions are arranged in a matter (or chiral) multiplet, 
vector bosons and fermions are arranged in a gauge (or vector) multiplet. 

Supersymmetry breaking terms, which do not reintroduce quadratic di- 
vergences, are added to the Lagrangian by hand |40j (in the notation of [H]): 

£ = -^(MaA'^A'^ + c.c.) - (m2);.</^*0, - (^U^<P4, + ^a*^V^0,0fc + c.c\ 

(2.3) 
There is a gaugino mass term (Ma) for each gauge group, scalar mass terms 
((?n^)*, U^) and a trilinear coupling a^^^ . These terms are added ad hoc 
to parametrize the current ignorance of the mechanism of supersymmetry 
breaking. 

If supersymmetry is gauged locally, gravity can be incorporated, leading 
to supergravity. The graviton, a massless spin-2 particle, must have a su- 
persymmetric partner (gravitino, spin-|). The gravitino can help solving the 
problem of generating the soft supersymmetry breaking masses flavor-blind 
where a hidden sector is introduced which communicates with the visible 
sector via gravitational interactions (gravitino condensation in the hidden 
sector). The hidden sector is connected to the gaugino and scalar soft su- 
persymmetry breaking terms. Additional terms ensure that different soft 
breaking masses can be generated for each superparticle. 



11 



2.2 The Minimal Supersymmetric Extension 
of the Standard Model 

Each Standard Model particle (or rather degree of freedom) receives a new 
particle as partner. The following section follows closely the notation of [TH 
The parameters of the MSSM are defined at the electroweak scale. 



spin-0 
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spin-1 


Qr, Ql 


q 
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g 


4,4 
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H, Z 
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h A 


h 




H± 


H±, W± 


W± 



Table 2.1: The (s)particle content of the MSSM 



Two Higgs doublets of opposite hypercharge are used to avoid triangle 
anomalies and FCNCs. Different Higgs fields give masses to the up- and 
down-type particles. The complete list of (s)particles and their superpartners 
is shown in Table [2?T1 



V = {ml + \fi\'^)Hi*Hi + {ml + \fi\'^)Hi*Hi-ml^{eijHiHi + h.c. 



Ti\2 



(2.4) 
Up-type fermion masses are proportional to fi, the vacuum expectation value 
of the first doublet, and down-type fermion masses to f2, the vacuum expec- 
tation value of the second doublet. The parameter /i is the supersymmetric 
Higgs mass parameter. The parameters governing the Higgs sector are usu- 
ally taken to be the mass of the lightest neutral CP-even Higgs boson (h) or 
CP-odd Higgs boson (A), tan/3 and mw from the electroweak sector. 

At tree level at least one Higgs boson should be lighter or as light as 
the Z boson. However, this stringent bound for the neutral Higgs bosons is 
loosened substantially by radiative corrections, which are driven mostly by 
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the heavy top quark mass and the top squark masses: 

, n 3q^ mf , ni2 mi 

Am^ = -^^ln^^ 2.5 

87r2 m^, m^^ 

The only strong prediction remaining after including these corrections is 
that the lightest neutral Higgs boson (h) must have a mass less than about 
140 GeV. 

In the leptonic sector each charged lepton receives two scalar partners, 
called L- and R-sleptons (Sl, 6^), where "L" and "R" refer to the chirality 
state of the corresponding lepton. For the neutrino there is one sneutrino 
(z/l)- In the gauge sector the partners of the neutral gauge and Higgs bosons 
mix to form the neutralinos: Xi)X2?X35X4- In turn the partners of the charged 
gauge and Higgs bosons mix to form the chargino mass eigenstates Xi ^ind 
xf- The gluon, as a massless spin-1 particle is grouped with the gluino, its 
fermionic partner. 

R-parite |12] is used to avoid certain yukawa terms {Xijk, X'ijk, Kjk) which 
are allowed in a general supersymmetric Lagrangian, but can lead to a rapid 
proton decay inconsistent with experimental limits. The conservation of this 
multiplicative quantum number in the following studies has as consequence 
that supersymmetric particles must be produced in pairs. They decay, possi- 
bly via cascade decays, until the lightest supersymmetric particle is reached, 
which is stable. 

The lightest neutralino is a good candidate for the lightest supersymmet- 
ric particle. It is color-neutral, not charged and interacts only weakly with 
matter. The lightest neutralino is a good candidate to solve the dark matter 
problem of the universe. 

Certain relationships among supersymmetric particles hold: 

m| = m? — cos^ 6'wni| cos 2/3 (2.6) 

Consequently the left-handed slepton must always be heavier than the sneu- 
trino, if tan fi is greater than unity. Similar relationships hold in the squark 
sector. 

The mass matrix for down-type squarks and charged sleptons can be 
written in the following way in the basis of left-handed and right-handed 
sfermion: 

m? mf ( Af — /i tan /3) 

; (2.7) 

mf (A? — u tan 13) m? 
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the off-diagonal matrix elements are proportional to the trilinear coupling Aj 
and to the fermion's mass. Since the mixing angle is governed by the fermion 
mass, it is expected that in the first two generations the weak eigenstates are 
equal to the mass eigenstates, so the trilinear mass terms for the first two 
generations are usually neglected in the scalar sector. For the up-type squarks 
tan/3 must be replaced by cot/3. 

The MSSM is governed by about 120 parameters. Not allowing CP viola- 
tion, requiring the trilinear terms to be real, requiring the soft supersymmetry 
breaking parameters in the scalar sector to be generation-diagonal to avoid 
additional sources of flavor changing neutral currents, reduces the number 
of parameters substantially. Additionally, the trilinear mass terms of the 
first two generations are neglected as their impact on the mass mixing is 
essentially undetectable experimentally. 

Another experimental requirement on the soft mass breaking terms in 
the first and second generation squark terms is derived from the inability 
to flavor tag, i.e., differentiate between u-, d-, s- and c-quarks. Adding 
everything together, a phenomenological MSSM is constructed. This pMSSM 
is governed by the following parameters in the supersymmetric sector and will 
be called MSSM in the following for simplicity: the Higgsino mass term (//), 
the mass of the CP-odd Higgs boson mass (uia), the ratio of the vacuum 
expectation values (tan/3), the gaugino mass parameters (Mi, M2, M3), 6 soft 
breaking parameters for sleptons (Mg^, Mgj,^, M^^^ M^Ar' ^n.' ^m)' ^ soft 
breaking parameters for the first and second generation squarks (Mq^^, Mqj^), 
3 soft breaking parameters for the third generation squarks (M^^^^, M^^ , Mqg^) 
and 3 trilinear parameters in the scalar fermion sector (A^, At, Ab). 

In addition to the supersymmetric parameters, for a precise prediction of 
supersymmetric observables, the following Standard Model parameters are 
necessary: the electromagnetic coupling constant at the Z-mass (a(mz)) the 
Fermi constant, the QCD coupling constant at the Z-mass (as), the Z-boson 
pole mass, mass of the bottom quark in the MSbar scheme, the pole mass 
of the top quark and the pole mass of the tau lepton. Thus in total 20 + 7 
parameters govern the MSSM. 

The electroweak scale where the supersymmetric parameters are defined 
needs to be specified. It is sometimes taken to be the geometric mean of the 
stop quark masses. This scale would then, in some scenarios, depend on the 
mass of a sparticle which has not yet been measured. Therefore the proposal 
of the SPA project [43] is followed by fixing the scale to 1 TeV/c^. 
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2.3 Minimal Supergravity 

While the MSSM is a general model depending on few assumptions, for 
practical purposes it is useful to define also a restricted model, motivated by 
supergravity: mSUGRA. The mSUGRA parameters are mostly defined at 
the GUT scale. The parameters are then evolved via renormalization group 
equations (RGEs) to the electroweak scale at 1 TeV. Every mSUGRA pa- 
rameter set has a corresponding MSSM equivalent, however not every MSSM 
parameter set has an mSUGRA equivalent. It is therefore interesting to study 
both approaches. 




Figure 2.1: Running of Oj, ; 
in the MSSM (full lines). 



2 4 6 8 10 12 14 16 18 
Log,o(Q/1 GeV) 

1, 3 in the Standard Model (dashed lines) and 



The Standard Model couplings of the three gauge groups SU(3), SU(2) 
and U(l), calculated according to the RGEs as a function of the mass scale, 
do not intersect in one point, as shown in Figure 1271) (a~^ = gl/{4:n)). In the 
MSSM ^7\ supersymmetric particles introduce an intermediate scale, modi- 
fying the slope after passing the threshold. In Figure \2A\ (taken from [Ti]). 
the three couplings intersect at at the GUT scale {Mu ^ 10^^ GeV). Inspired 
by flavor neutral gravity, it is assumed that the gaugino mass parameters are 
also equal at the GUT scale. As a consequence Mi, M2 and M3 are reduced 
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to one parameter at the GUT scale: nii/2- The following relationship holds 
at any scale: 

Ml _ M2 _ Ms _ mi^ , . 

91 92 9i 9v 

gu is the coupling at the unification scale, gl is related to the strong cou- 
pling constant, and the couplings gi, g2 are related to the usual electroweak 
constants hj g2 = g and gi = ^/5/3g' 

Additionally one assumes that the soft supersymmetry breaking terms for 
the sfermions are also universal at the GUT scale (mo). This, however, does 
not imply that all sfermions have equal mass at the electroweak scale. Ra- 
diative corrections for left-handed sfermions are different from right-handed 
sfermions, since their couplings to gauge bosons are different. The solutions 
of the renormalization group equations for the masses of the right-handed 
slepton, the left-handed slepton and the sneutrino at the electroweak scale 
take on the following form: 

m? = mg + 0.22M|-sin2^wm|cos2/3 + m^ 

m? = m2 + 0.75M2-0.5(l-2sin2^w)m|cos2/3 + m2 (2.9) 

m? = m2 + 0.75M2 + 0.5m|cos2/3 

A universal trilinear coupling is assumed at the GUT scale and one re- 
quires that the electroweak symmetry breaking is driven by supersymme- 
try [H]. This determines the supersymmetric Higgs parameter fi up to its 
sign. 

Thus in the supersymmetric sector the following parameters remain: the 
ratio of the Higgs vacuum expectation values defined at the electroweak scale 
(tan/3), the sign of fi in the Higgs sector, the universal gaugino breaking 
parameter (mi/2), the universal sfermion breaking parameter (uiq) and the 
trilinear term (Aq). The last three parameters are defined at the GUT scale. 
Additionally, as before for the MSSM, the Standard Model parameters have 
to be specified. 



2.4 Benchmark points 

Given the number of parameters, supersymmetry can lead to a wealth of 
different signatures at colliders. To enable comparisons between experiments 
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Table 2.2: Summary of mSUGRA parameter sets with similar collider phe- 
nomenology, /i is positive for all cases. 



and/or colliders, sets of parameters are defined to represent typical signa- 
tures. Most commonly known is the point Snowmass Points and Slopes 
(SPS) la [45j which has been studied in detail in the LHC-ILC report [16]. 
Points with similar phenomenology have been studied in ATLAS (SU3) and 
CMS (LMl [47J). A summary of the points is shown in Table 12.41 Com- 
mon to all points are relatively light supersymmetric particles, thus an early 
discovery at the LHC would be possible. 
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Figure 2.2: The particle spectrum of 
the SPSla point is shown. 

In the following SPSla will be studied as an example. The resulting 
spectrum of (s)particles is shown in Figure 12.21 The parameters of this 
point lead to first and second generation squarks with masses of the order 
of 500 GeV/c^. The gluino has a mass of about 600 GeV/c^. The lightest 
right-handed sleptons have a mass of 150 GeV/c^, their left-handed counter 
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parts have masses of 200 GeV /c^. The hghtest Higgs boson mass is at the 
edge of the LEP bound. There is a group of hght neutrahnos and chargino 
at about 200 GeV/c^ and below as well as a group of heavy chargino and 
neutralinos. The right-handed squarks decay predominantly directly to their 
standard model partner and the lightest neutralino. The left-handed squarks 
have a long decay chain via the second lightest neutralino and a right-handed 
slepton open. As the lightest stau (with a left-handed component) is lighter 
than the second-lightest neutralino and the lightest chargino, the chargino 
decays with a branching ratio of almost 100% to stau and neutrino. The 
neutralino decay to stau and tau has a branching ratio of 90%. 



Chapter 3 

Colliders and Experiments 



The Large Hadron Collider (LHC) at CERN has started operations on Septem- 
ber 10, 2008. The LHC will collide proton beams with a center-of-mass en- 
ergy of up to 14 TeV. Two multi-purpose detectors, ATLAS and CMS, are 
located on the LHC ring. The ILC, an e"'"e~ collider with a center-of-mass 
energy of 500 GeV, extendable to 1 TeV is awaiting approval pending the 
results of the LHC. The ILC is expected to start operations eight years after 
approval. 

In this chapter the LHC is briefly described. The ATLAS detector is 
discussed with emphasis on its performance, especially for electron and pho- 
ton reconstruction. A hardware aspect of the readout of the Liquid Argon 
calorimeters is analyzed in more detail. The ILC and an associated detector 
concept is discussed last. 

3.1 LHC 

The LHC is described in detail in Ref. j48j. The nominal mode of operation 
of the LHC is to provide a luminosity of lO^'^cm'^s"^ at a center-of-mass 
energy of 14 TeV colliding proton beams. The bunch crossing frequency will 
be 40 MHz. The beams circulate in a tunnel with a circumference of 26.7 km 
which is connected to the CERN accelerator complex via two 2.5 km long 
transfer lines. 

The luminosity at a collision point is described by the following formula: 

^^iVX^evT^ (3.1) 
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where 7 is the relativistic gamma factor (for protons at 7 TeV about 7000), e 
the normahzed transverse beam emittance (3.75/im), /3* the beta function at 
the interaction point (0.5 m), /rev the revolution frequency, Nb is the number 
of particles per bunch (maximum 1.1-10^^) and rib the number of bunches per 
beam, nominally 2808. Note that Nb enters as a square while Ub only linearly 
as the luminosity is defined per interaction point. Bunches are separated in 
the beam pipe by 7.5 m. To minimize beam collisions in the 130 m common 
beam pipe in the interaction regions (corresponding to about 17 unwanted 
collisions for each bunch of a beam), a small crossing of 150-200 /irad will 
be used. F is the geometrical reduction factor due to the non-zero cross 
angle. The size of a bunch is expected to be about 15 /im in the direction 
transverse to the beam axis and about 7 mm in the direction of the beam 
axis (z-direction). 

The LHC tunnel consists of eight straight sections and eight arcs with 
a smaller radius, the LHC magnets are built to reach a strength of up to 
9 T with 8.33 T being the nominal operating field for a beam momentum of 
7 TeV/c. A total of 1232 superconducting dipoles have been built in industry. 
Adding in quadrupole and other beam equipment, the LHC consists of about 
4000 elements. 

The dipoles are operated at a temperature of 1.9 K. This low temperature 
is achieved in a two stage process, first by cooling down to about 90 K with 
liquid Nitrogen and then cooling down to the operating temperature with 
super-fluid liquid Helium. The LHC has been divided into eight sectors. The 
cool-down of one sector to operating temperature takes about two months. 
The energy that will be stored in the beams will be around 362 MJ (60 kg 
of TNT) and an additional 600 MJ is stored in the magnet system. A fast 
system is able to extract the beam if necessary and dilute it before reaching 
an absorber. 

The injection sequence of the beams is shown in Figure [3?T1 The protons 
are first accelerated to 50 MeV/c in a linear accelerator. In the Proton Syn- 
chrotron Booster (circular) they reach 1.4 GeV/c before being transferred 
to the Proton Synchrotron (PS). In the PS the protons are accelerated to 
25 GeV/c and injected into the Super Proton Synchrotron (SPS). The pro- 
tons with an energy of 450 GeV are then transferred from the SPS to the 
LHC. The beam is captured and accelerated by a superconducting RF sys- 
tem. In contrast to electron beams, the synchrotron radiation loss for protons 
is small ~ (me/nip)^. The energy gain per turn during the ramp to 7 TeV is 
485 keV. The expected beam lifetime (until the luminosity has decreased to 
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Figure 3.1: The CERN accelerator complex from the Linear accelerator to 
the LHC. 



1/e) is expected to be about 15 h. The average filling time should on average 
be about 7 h. 

The integrated luminosity per year is expected to be 10 fb^^ running 
at lO^^m^^s^^. For the nominal luminosity of 10^^m~^s~^, an integrated 
luminosity of 100 fb~ is expected. 

Beams were injected in both rings on September 10, 2008 and circulated 
many times. The RF-system has captured the beams leading to a nearly 
infinite theoretical lifetime. However, during the commissioning of a sector 
to the highest beam energy, a superconducting soldering joint turned normal 
conducting with a resistance of a few nQ. The small resistance was large 
enough to heat the super-fiuid liquid helium turning it into its gaseous state. 
The volume expansion was so quick that the safety valves did not follow 
quickly enough. The shock wave displaced several magnets. About 50 ele- 
ments of the machine have been repaired or replaced by spares. Additional 
safeguards are being put in place. 

After the incident in 2008, repairs of the cryogenic as well as the magnet 
system are underway. The current startup plan aims to restart the machine 
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in November 2009 with a maximuiii center-of-mass energy of 7 TeV. The 
ramp up to the nominal center-of-mass energy of 14 TeV is foreseen for 
later. 



3.2 The ATLAS Detector 
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Figure 3.2: View of the ATLAS detector. 



The ATLAS detector |19], as shown in Figure [3^ is a typical Air detector. 
Its innermost parts consist of the tracking devices in a 2 T solenoidal magnetic 
field covering the pseudo-rapidity [t] = — logtan6'/2, 6 is the polar angle 
with respect to the beam axis) of \ri\ < 2.5. Electromagnetic and hadronic 
calorimeters follow further away from the interaction point covering rj down 
to 4.9. The muon detectors are further out in toroidal magnetic field covering 
|?7| up to 2.4. 

At the LHC the total inelastic cross section for proton-proton interactions 
is 80 mb. At nominal luminosity every event will be accompanied by an 
average of 23 minimum bias events. The LHC will collide protons every 
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25 ns in the center of ATLAS. The high cross section necessitates a fast 
readout of the detector as well as radiation hard or tolerant electronics for 
components on the detector. 

The tracking detectors are immersed in a solenoidal field of 2 T. It is 
expected that 1000 particles will come from the collision point with a fre- 
quency of 40 MHz within a pseudo-rapidity of 2.5, thus the track density 
will be large. The innermost tracking device is the Pixel detector, which is 
segmented in R-0 and z. Each track crosses three pixel layers, the innermost 
is at a radius of 5 cm. The intrinsic accuracy in the barrel is 10 /im in R-0 
and 115 /zm in z from 80.4 million readout channels. 

The Pixel detector is followed by the SCT, consisting of silicon strips with 
a slightly coarser granularity. The intrinsic accuracy in the barrel is 17 /im 
in R-0 and 580 fim in z. There are 6.3 million readout channels for the SCT. 

The TRT provides more points for tracking than the Pixel and SCT, 
however each point has less precision. On average 36 points are provided for 
tracking with R-0 information with an intrinsic accuracy of 130/im per straw 
up to |?7| < 2. Yet no information is available in the z direction. The TRT 
has 351000 channels. In addition to its tracking capabiities the TRT also 
provides particle identification. The TRT measures the transition radiation 
which is expected to yield a higher signal for electrons than for pions. This 
information is used for the identification of electrons described later in this 
Chapter. 

The resolution for the tracking detectors has been determined in a re- 
alistic simulation of the ATLAS detector. Using the parametrization ax = 
crx(inf)(l©px/PT), where crx(inf) is the asymptotic resolution at infinite mo- 
mentum and px is a constant for which the intrinsic and multiple-scattering 
terms are equal. The resolution in the barrel in simulation correspond to 
cTx (inf) = 0.34 TeV~^ and px = 4:4: GeV. For a transverse momentum of 
100 GeV/c, the relative transverse momentum resolution pt-(j{1/pt) is about 
0.04. 

The electromagnetic calorimeter of ATLAS is a sampling calorimeter. 
The passive material is lead (and stainless steel for the mechanical stability). 
The active material is Liquid Argon, operating at a temperature of about 
90 K. The particularity of the design is the accordion structure as shown 
in Figure 13.31 ensuring a coverage in without dead zones. The depth of 
the calorimeter varies roughly, as function of r/, from 25 Xq to 30 Xq. An 
electrical field of 2 kV is applied in the 2 mm Liquid Argon gap. 

The calorimeter has three segments in depth providing rj coverage up to 
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Figure 3.3: View of a module of the accordion calorimeter of the ATLAS 
barrel during assembly. 



3.2 in addition to a presampler (PS) covering \ri\ < 1.7. The first segment, 
the strips, typically has a cell size of A?^ = 0.003, coarser in (p with a size 
of 0.1. The fine granularity allows to separate efficiently photons from vr". 
The main part of an electromagnetic shower is deposited in middle section 
where the typical readout cell size in t] and direction is 0.025 x 0.025. The 
thin back section has the same cell size in as the middle and is coarser in 
1] with a cell size of 0.05. The expected energy resolution for electrons and 
photons is 10%/-\AE and a global constant term of 0.7 % must be achieved. 
The linearity is required to be of the order of a per mil or better to measure 
precisely the absolute masses of particles. 

The hadronic and forward calorimeters in the endcap are also based on 
Liquid Argon technology. The absorber of the hadronic calorimeter is made 
of copper. The calorimeter has four segments in depth. Their size in Ar/ x Acp 
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is 0.1 X 0.1 in the region of |?7| < 2.5 and 0.1 x 0.2 at larger t] for a total 
of 5632 readout channels. The forward calorimeter (FCAL) consists of three 
high-density modules. The first one, for electromagnetic measurements, is 
made of copper, the other two are made of tungsten. Concentric rods and 
tubes are installed parallel to the beam axis. The gap between rod and tube, 
filled with Liquid Argon as the sensitive material can be as small as 0.25 mm. 
This ensures a fast signal in the forward region where the pileup effects are 
expected to be stronger. 

The tile calorimeter is situated behind the electromagnetic calorimeter in 
the barrel. Scintillating tiles provide the signal of the sampling calorimeter 
while steel is used as absorber. Three layers are provided in depth. The 
readout cells are of size 0.1 x 0.1 in the first two layers and 0.1 x 0.2 in the 
last layer. Measured in interaction lengths, the tile calorimeter is almost lOA 
deep at ?7 = 0. 

The muon system provides coverage for |?7| < 2.7. The muons are de- 
flected in large air core toroids. 1800 hall sensor monitor the magnetic field. 
Monitored drift tubes are used over most of the r] range to provide infor- 
mation for the track reconstruction, cathode strip chambers are used in the 
endcap due to the higher occupancy. Resistive plate chambers (barrel) and 
thin gap chambers (endcap) are used for triggering. 

ATLAS uses a three stage trigger system to reduce the rate from 40 MHz 
in level-1 to 75 kHz and further down to 200 Hz in the two high level trigger 
systems using the full detector information. 

3.2.1 Electromagnetic Test beam 

A vigorous test beam program has been pursued by ATLAS. Components 
of all sub-detectors have been tested. A pre-series module of the electro- 
magnetic calorimeter was tested in the test beam with electron beams of 
fixed energy from 10 GeV to 245 GeV. Four barrel modules, correspond- 
ing to 12.5% of the barrel calorimeter and three endcap modules (19% of 
the endcap calorimeters) which are now installed in ATLAS have also been 
tested. 

The main results of this test beam campaign of nearly five years (1998- 
2002) have been reported in Refs. [501 EH [53 [531 |Ml EH ES] . The uniformity, 
i.e., the variation of the calorimeter response for a fixed energy of 245 GeV, 
in a barrel module of effective size Arj x A(f) = 1.2 x 0.2 is better than 0.5%. 
The uniformity of the calorimeter is especially important for the measure- 

25 



ment of thin resonances, e.g., the Higgs boson decaying to two photons, or 
sharp endpoints (the lepton-lepton invariant mass in a supersymmetric de- 
cay chain). If the intrinsic width is smaller or about equal to the uniformity, 
the width of the distribution, e.g., of the invariant photon-photon mass, will 
be dominated by the detector effects. The result obtained is in excellent 
agreement with the requirements set forth before the construction of ATLAS 
started: uniformity in a region of 0.2 x 0.4 better than 0.5%. The number 
shows the excellent control of the long and difficult production process of all 
components of the electromagnetic calorimeter. 

In addition to the uniformity the energy resolution with its sampling term 
and local constant term contributes to the width of the measured invariant 
masses and edges. In the test beam the energy resolution was measured to 
be (at a fixed r] position) better than 10%/vE and the local constant term 
to be better than 0.3%. The performance expected for the electromagnetic 
calorimeter of ATLAS has been achieved. 

While uniformity and energy resolution determine the width, the energy 
scale (or central value) needs to be known precisely. Therefore the linearity 
of the calorimeter is very important. The linearity was measured to be of 
the order of 2 per mil in a range from 20 GeV to 180 GeV in a special test 
beam setup. 

The test beam has been tremendously important to develop the methods 
for reconstruction of electrons and photons. However the results obtained 
with the modules now installed in the ATLAS detector cannot be transferred 
directly from the test beam to ATLAS, because the material distribution in 
the test beam is quite different from the one in situ in ATLAS. The strategy 
chosen was to develop the methods on the test beam, optimize the test beam 
Monte Carlo, heavily used for the linearity study, and apply the methods to 
the simulation of the full ATLAS detector. The Monte Carlo provides the 
(indirect) link between the test beam and the (future) ATLAS data. 

In order to have a reliable measurement from the electromagnetic calorime- 
ter it is also necessary to monitor every aspect (stability of the electronics 
system etc). The dependence of the energy response as function of the liquid 
Argon temperature has been measured to be 2%/K for a fixed beam energy 
of 245 GeV. The temperature stability in the test beam cryostat was shown 
to be better than 7 mK. The Front End electronics, in particular the Front 
End boards (FEBs) where the signal is amplified, formed and digitized, were 
tested and their design was validated. 
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3.2.2 The Front End Board of the LAr Calorimeters 

The Front End Electronics of ATLAS consist of several components. Here, 
as an example the Front End Board (FEB), described in detail in Ref. |33] is 
discussed. The FEB is located at large radius on the cryostat in the transition 
region between barrel and endcap. It amplifies the signal of the Liquid Argon 
detectors of ATLAS, with the exception of the HEC, where the amplifiers 
are located inside the cryostat. 

After preamplification the FEB forms the triangular ionization signal into 
a bipolar shaper where the peak of the signal is reached typically after 50 ns. 
This corresponds to only 2 bunch crossings, whereas the ionization signal 
has a duration of approximately 450 ns. Additionally the bipolar shaping 
allows to treat the minimum bias/pileup events, which are typically at low 
transverse momentum, as an additional noise contribution to the resolution 
of the electromagnetic calorimeter. The shaper also serves as a first step 
summing device for the trigger. The analog signal is then sampled at 40 MHz 
and saved in a switch capacitor array (SCA) for about 1.5 /is until the first 
level trigger decision is received. Upon a trigger, the signal is digitized and 
sent to the Read Out Device (ROD) via an optical link. 

The FEB treats 128 readout channels of the calorimeter. 1627 FEBs were 
produced and configured in 19 different flavors. The production, assembly 
and test was shared among three labs: NEVIS (University of Columbia), 
Brookhaven National Laboratory (BNL) and LAL. While NEVIS was re- 
sponsible for the production and the first digital tests, LAL and BNL each 
configured, tested and assembled 50% of the FEBs. The operation lasted 
for almost a year. Three engineers, three technicians and a physicist were 
involved in the development, maintenance, improvement and functioning of 
the test setup up at LAL alone. This allocation does not cover the 10 years 
of R&D of the FEB components in radiation tolerant technology (see [33] for 
more details). 

At LAL, the FEBs were tested one by one in specially developed test 
stand. The shaper output is in three gains. First the pedestals and the 
electronic noise were measured. The distribution of the noise RMS of the 
high gain channel is shown in Figure l37i] for all FEBs of a flavor. This test 
allowed to detect dead channels. 

To test the signal integrity, a calibration pulse was sent through the elec- 
tronics chain and the signal was recorded in steps of 1 ns with interleaved 
events as shown in in Figure 13.51 The analysis of the signal form yielded in 
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Figure 3.5: Reconstructed signal of 
the FEB in medium gain as function 
of the time (from [33] ) . 



particular the time of the signal peak. The peak was required to fall into an 
acceptance window 7 ns long. The minimum and maximum of the time win- 
dow were fixed to be the same for all FEBs of the same preamplifier flavor. 
A wrong time constant of the shaper, usually observed as a signal reaching 
its maximal amplitude before the minimum time required, was detected this 
way. 

The distribution of the measured linear part of the gain is shown in Fig- 
ure 13.61 The slope is sensitive for the detection of dead channels, accelerated 
shaper timing and channels with linearity problems. In normal running in 
ATLAS only one of the three gains is read out. The gain choice is made 
by a dedicated chip (GSEL) on the FEB. The integrity of the choice, cor- 
responding in practice to the setting of two thresholds in the medium, was 
verified. 

The FEB provides the signals for the trigger level-1. In order to avoid 
excessive triggering due to noisy channels, individual calorimeter readout 
cells can be turned off. Figure 13.71 shows the test for one FEB as example. 
This test also checks the correct signal routing of the trigger signals on the 
FEB. Additionally the gains as well as the shapes of the trigger outputs were 
measured and compared to references as well. 

On the 786 boards configured and tested at LAL with roughly 75000 
shapers, SCAs and Preamplifiers, 801 repairs had to be performed. Part of 
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Figure 3.7: Test of the FEB trigger 
output. No signal is expected for the 
first 32 entries and for the next 128 
entries a large signal should be ob- 
served (from [33]). 



the repairs (faulty chip replacement) was performed in industry, another part 
in house. 

The FEBs are cooled with water circulating in plates which are mounted 
individually on each FEB. To test that the system has no leaks, an air pres- 
sure test at 3 bar was performed. The water system in the experiment is an 
under-pressure system which ensures that, in case of a leak, air would enter 
the system, but no water would flow out. The FEBs were shipped to CERN 
after a further functional verification. 

Since the completion of this operations the test bench has been kept 
functional. It was used again in 2008/2009 for a small production of about 
40 additional boards. After more than a year of operation on the detec- 
tor several boards were malfunctioning due to corrosion. The corrosion was 
traced to an error in the cleaning procedure after chip replacement at the 
subcontractor. In addition the peaking time of the signal was unstable over 
time. The FEBs were removed from the detector, cleaned if necessary, or 
replaced. The peaking time of the shaper was stabilized by cutting pin con- 
nections between the shaper and FEB board (1627 boards times 32 circuits 
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per board). 

The production, assembly and test of the FEBs was a difficult endeavor. 
It was a successful operation as the time delay given by the collaboration 
was respected. A large number of problems on the boards were detected and 
repaired before installation. 

3.3 Performance of the ATLAS Detector 

The ffist physics results expected from collision data will be essential to 
prepare the physics calibration and measure the performance in situ. In 
order to evaluate the current understanding of these issues, physics studies 
and a test of the analysis model, ATLAS has performed a Computing System 
Commissioning (CSC) evaluation [31] . 

For the CSC studies more than 25 million events were simulated, digitized 
und reconstructed. The event samples comprised single particles, essentially 
for calibration and identification, as well as complex physics events at the 
nominal center-of-mass energy of 14 TeV. The best known description of the 
detector, including the experience gained from the test beam studies of all 
sub-detectors, was used in the simulation. 

Muons are reconstructed with an efficiency of more than 95% for tracks 
with transverse momenta of more than 5 GeV/c. Good coverage in rj is 
obtained in combination with a low fake rate as shown in Figure ESI (Left). 
Muons can be reconstructed in the muon spectrometer alone as well as in 
the inner detector. The combination of the two significantly reduces the fake 
rate by an order of magnitude. In order not to rely on Monte Carlo for the 
determination of the efficiency, the production and decay of the Z boson will 
be used. With an integrated efficiency of 100 pb~^, using the well known Z 
boson mass, a statistical precision of 0.1% and a systematic error of 1% can 
be obtained. With only half as much integrated luminosity the muon energy 
scale will be known with an accuracy of ±0.5 GeV for transverse momenta 
of 50 GeV/c. 

Jet reconstruction is performed in several stages. First clusters are re- 
constructed in the calorimeters and calibrated locally [21]. Jets are then 
reconstructed with a jet algorithm, most frequently the cone and the kt algo- 
rithms though others are also under study. Noise, pileup and the underlying 
event are corrected for as the last step. The jet energy scale, as shown in 
Figure [XH is correct at the percent level. Photon plus jet events and Z boson 
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Figure 3.8: (The efficiency for recon- 
structing muons is shown as func- 
tion of the transverse momentum for 
the muon spectrometer reconstruc- 
tion, combined reconstruction and 
combined with the use of segment 
tags (taken from [l9]). 



Figure 3.9: The signal hnearity of 
jets reconstructed with Rcone = 0.7, 
in two different regions of \ri\ is 
shown as a function of Etruth for 
jets at the electromagnetic energy 
scale and fully calibrated jets (taken 
from [19]). 



plus jets will be used to carry out the in situ calibration. Due to the large 
production cross sections, only 10 pb~^ will be necessary to reach a percent 
precision for low (less than 80 GeV) transverse energies. 

In the E'™** reconstruction the energy deposits in the readout cells of 
the calorimeters are refined by the identification of muons, electrons and 
jets. Not using this knowledge would lead to an offset of about 30% with 
respect to the true missing energy. The linearity is about 5% and the EJp^^^ 
resolution is expected to be 0.57- ^/SEt, as shown in Figure 137701 The scale 
will be verified for low EJp^'^^ with the Z boson decaying to taus. 

The hadronically decaying taus are either reconstructed by a tracking 
based or a calorimeter based algorithm with an efficiency of 30% for a jet 
rejection of 10'^. Clean samples of taus will come from the production of the 
W and Z bosons. The visible mass reconstructed from a leptonically and a 
hadronically decaying tau will peak at about 54 GeV/c^. This can then be 
used for the calibration of the hadronic r decay. 

For b-tagging a working point at 60% efficiency is defined, where the 
rejection of light jets is about 30 for the simplest algorithm. Using the 
longitudinal and transverse impact parameter significance, the rejection is 
increased by a factor two. Reconstruction of the secondary vertex leads to 
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Figure 3.10: The resolution of the 
j^miss yg(.^Qj, jg shown for QCD di-jet 

events and Higgs bosons (A) includ- 
ing a parametrization 0.57i/SEx. 
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Figure 3.11: The expected differen- 
tial spectrum for single jets is shown 
after the level-1 trigger and pre-scale 
factors as well as without (taken 
from 



a light jet rejection of 150. The b-tagging performance will be calibrated in 
situ with top quark pair production to a precision of 5-10% and with muon 
plus jets at low momentum. 

Trigger menus have been developed for different values of the instan- 
taneous luminosity, expected to be several orders of magnitude below the 
nominal 10^'^cm~^s~^ at startup. The total rates are within budget and the 
turn-on curves are reasonably sharp. Differential jet production will be mea- 
sured over many orders of magnitude, as shown in Figure [3.111 

3.3.1 Electron and Photon Reconstruction 

The reconstruction, calibration and identification of electrons and photons 
are an example of the transfer of test beam knowledge to ATLAS simulation 
and reconstruction. 

In the standard algorithm for isolated electrons and photons, clusters 
are reconstructed from calibrated readout cells. The cluster algorithm uses a 
fixed size, independent of the energy of the cluster. A track is then matched to 
the cluster geometrically and uses the ratio of calorimetric energy to tracker 
momentum. To qualify the candidate as electron candidate the track is 
required not to be associated to the reconstructed conversions. Conversions 
are reconstructed independently in the tracking system alone. Clusters with 
conversion tracks and/or clusters without associated tracks are classified as 
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photon candidates. 

For the cahbration the difficulty hes in the simultaneous optimization of 
resolution and linearity (energy scale). Two methods of electron and photon 
calibration have been developed to address this problem. The ffist one is 
a weighting technique which uses the presampler and the energy deposited 
in the calorimeter segments to derive from Monte Carlo the correction fac- 
tors to be applied. These factors correct for the energy loss upstream of the 
calorimeter. The material in front of the calorimeter (inner detector, services, 
mechanical structures) is several radiation lengths, leading to bremsstrahlung 
and photon conversions. The second method uses the Monte Carlo to cor- 
relate the energy lost with the shower depth and derive a correction (dead 
matter calibration). The application of these methods leads to a linearity 
of the order of per mil up to several hundred GeV. In addition to this lon- 
gitudinal weighting, local corrections are also applied to correct the cluster 
position as well as to correct for the local energy variations. For details see 
Ref. |3T]. 

As the cross-sections for QCD jet production are several orders of mag- 
nitude larger than, for example, the production of the Z boson decaying 
to electrons, a high jet rejection is necessary while maintaining high effi- 
ciency. The electron and photon identification is based on the use of the 
shower shapes, calculated from the order 100 readout cells building an elec- 
tromagnetic cluster, from tracking information and from track-matching in- 
formation. In particular on the tracking side the TRT transition radiation 
measurement is used to differentiate between charged pions and electrons. 
On the calorimeter side the fine granularity of the ffist calorimeter segment 
(strips) allows to reduce the pion (decaying to two photons) by a factor three. 

The differential jet cross sections drop sharply as function of the trans- 
verse energy. The cross section for some discovery channels of new physics 
might be large. On the other hand, sensitivity to rare new physics pro- 
cesses with electrons and photons must be retained. These different needs 
led ATLAS to define three electron identification qualities, one (tight) with 
the highest possible jet rejection of 10^, a second one (medium) with a re- 
jection of 2200 (with higher efficiency) and third one with a jet rejection of 
570. The corresponding efficiencies are shown in Figure 13.121 as function of 
the transverse energy. For photons a jet rejection of 8000 has been obtained 
for an efficiency of 85%. 

The Z boson data will be essential to determine the electron reconstruc- 
tion and identification efficiency in situ. One electron is well identified (tag) 
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Figure 3.12: The electron efficiency 
is shown as function of the trans- 
verse energy (taken from [31j) for su- 
persynimetric events as well as single 
particles. 
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Figure 3.13: The reconstructed 
Higgs mass is shown. The shaded 
histogram are the events where 
at least one photon has converted 
(taken from |31j). 



and the other electron is used to study the efficiency (probe). Good agree- 
ment is obtained between the tag&probe method and Monte Carlo truth. A 
statistical error on the efficiency of ±0.1% and a systematic error of ±1.5% 
can be obtained with an integrated luminosity of 100 pb~^. 

The Monte Carlo can only help to determine the best possible reconstruc- 
tion for the linearity. The absolute scale must be determined by using the 
process Z— )■ ee. The ATLAS electromagnetic calorimeter has a uniformity 
of 0.5% in regions of Ar/xA^ = 0.2x0.4, as shown in the test beam studies. 
These regions are inter-calibrated to obtain a long range constant term of 
about 0.4% with an integrated luminosity of 100 pb~^ Thus the goal of a 
global constant term of 0.7% is within reach. 

As an illustration of the performance of the electron and photon recon- 
struction the reconstructed invariant mass of the photons from the decay of 
the Higgs boson to two photons is shown in Figure 13.131 The energy scale is 
correct to about 500 MeV. The shaded area are events where at least one of 
the photons has converted. 

In addition to the standard algorithm described here, a second algorithm, 
more appropriate for the reconstruction of non isolated electrons, e.g., in the 
vicinity of jets, has also been developed by ATLAS. As the ffist source of 
electrons at the LHC will be the copious production of electrons from J/i/j/T 
decays with rather low transverse energy, it will be useful to use the two 
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algorithms to study the overlap and possible improvements. 

The lessons learned from test beam measurements have been transferred 
successfully to the reconstruction of physics events simulation with a detailed 
description of the ATLAS detector. The performance of the reconstruction, 
e.g., electrons and photons, is in agreement with expectations providing a 
good basis for future studies and most importantly the collision data from 
2009 on. 

3.4 e"'"e~ Linear Collider 

Two designs for a linear collider are being pursued currently: ILC and 
CLIC The ILC is an e"'"e~ collider design with a center-of-mass energy of 
500 GeV upgradable to 1 TeV [57]. While the nominal center-of-mass energy 
is lower than that of the LHC, all the energy is available in the collisions in 
contrast to the LHC. The properties of an e^e~ collider complement the 
LHC studies. Therefore the ILC is considered to be the next project for high 
energy physics, depending on the results obtained by the LHC. 

The layout of the ILC is shown in Figure 13.141 The tunnel length is 
about 31 km. Electrons and positrons are accelerated in tunnels of 11 km 
length. The beam delivery system is about 4.5 km long. The crossing angle 
is planned to be 14 mrad. Two detectors will be installed in a push-pull 
system, alternating their data-taking periods. 

The main technological challenge of the ILC is the acceleration of the 
electrons, whereas at the LHC the challenge is the magnet system. The ILC 
design calls for the use of superconducting cavities, operating at a frequency 
of 1.3 GHz, with an average gradient of 31.5 MeV/m. In addition to this 
challenge, it is planned to have an electron polarization of more than 80%. 
The positrons are produced by conversion of photons (passing the electron 
beam through an undulator) into e"'"e~ pairs. The current design of the 
machine has a positron polarization of 30%, which in a upgrade could be 
increased to 60%. The polarisation must be measured with a precision of 
0.1% via the compton effect. 

The center-of-mass energy can be varied from 200 GeV to 500 GeV. 
Short periods at 91 GeV (Z-resonance) for calibration will also be possible. 
The precision of the knowledge of the beam energy will feed directly into the 
absolute uncertainty of the measurement of any particle mass. It is expected 
to be able to measure the beam energy to 200 ppm by using two different 
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Figure 3.14: Layout of the ILC 
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methods. Upstream the deflection of the beam through a dipole field will 
be measured. Downstream the synchrotron radiation of the beam passing 
through a string of dipole magnets will be measured. The beam has a pulse 
rate of 5 Hz with a pulse length of about 1 ms. In each pulse there will 
be 2625 bunches with 2 ■ 10^° particles per bunch. The design luminosity of 
the ILC is 10^^m~^s~^. The RMS of the beam size at the interaction point 
is expected to be about 6 nm (vertical) and about 640 nm (horizontal). In 
addition to initial state radiation, 2.4% (RMS) of the energy will be lost due 
to beamstrahlung. The downstream beam energy measurement will mea- 
sure the beam energy spread as well as the beam energy distribution after 
beamstrahlung has occurred. 

While the preparation of the ILC is well advanced, R&D is also performed 
for a collider further in the future . A second generation e"'"e~ collider called 
CLIC (|60] and references therein) is being developed. The CLIC concept is 
based on a novel scheme decelerating one beam to accelerate another one. 
CLIC has the potential to reach a center-of-mass energy of 3 TeV. 




Figure 3.15: Event display of a simulated event in a detector for the ILC. 
The fine granularity of the calorimeter used in the particle flow is clearly 
seen. 

Currently there are three detector concepts under study, competing for 
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the two ILC detector spots. The two main ones are ILD (International Linear 
Detector) [SS] and SiD (Sihcon Detector) [SH]. They differ in the technology 
chosen for the tracking sub-detectors. SiD has silicon based tracking with a 
strong magnetic field of 5 T, thus choosing to build a compact detector. ILD 
has chosen a larger volume design with a lower magnetic field of 4 T and with 
a large gaseous TPC in addition to the silicon vertex detector. The goal for 
the tracking system is to reach a precision of o"(1/pt) = 10~^GeV~^ above 
10 GeV/c, an order of magnitude more precise than the precision achieved 
by the LEP detectors. 

The calorimeters chosen by both collaborations are similar: the electro- 
magnetic calorimeter is a very high granularity detector. For this reason it 
is sometimes called tracker-calorimeter, with 100 million readout channels 
(the ATLAS electromagnetic calorimeter has about 150000 readout chan- 
nels). The energy resolution for electrons is about 15%/a/E for electrons, 
which is less ambitious than LHC experiments (e.g. ATLAS 10% or CMS 
2%). However the main physics processes expected at the ILC are multi- 
jet processes where the individual performance is less important than the 
combined reconstruction of all sub-detectors in the particle flow concept. A 
jet resolution using Particle Flow of 30%/a/E is expected, improving on the 
50% expected in ATLAS. An event display showing the fine granularity of 
the electromagnetic calorimeter is shown in Figure 13.151 
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Chapter 4 

Determination of 
Supersymmetric Parameters 
and Higgs Couplings 



The determination of the fundamental parameters, be it supersymmetric pa- 
rameters or the Higgs couphngs, critically depends on detailed experimental 
simulations of measurements and errors at the LHC and at the ILC While 
the techniques described here are general, the well-understood parameter 
point SPSla [15] is used as a specific example. This point has a favorable 
phenomenology for both LHC and ILC. 

SPSla was invented as a parameter set for the study of the LHC and 
ILC discovery and measurement potential. The Standard Model electroweak 
measurements, the b-physics precision observables, [g — 2)^ [78] with the 
addition of the WMAP [77| measurement of the relic density already allows 
to delimit interesting regions of parameter space without the direct obser- 
vation of supersymmetric particles. In particular, as shown in [79l IHD] , the 
supersymmetric mSUGRA fit of these observables results in a preferred value 
for the lightest Higgs boson mass of mh = 113.3GeV/c^, pushing the Higgs 
boson mass closer to the limit of direct searches at LEP of 114.4 GeV/c^ [T8]. 
The effect of using the LEP exclusion is about 3 GeV/c^. The best-fit point 
is close, albeit with rather large errors, to the definition of SPSla [81]. This 
is a further motivation to study the phenomenology of this parameter set in 
detail. Using a different technique, the MSSM parameters have been studied 
using the precision observable as shown in Ref. [82] . 

Discovering supersymmetry, and in particular SPSla, at the LHC is con- 
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sidered to be feasible for masses up to 2.5 TeV/c^ (squarks and gluinos). In 
general, supersymmetry provides signatures with several high-pT jets and 
large missing energy. The search for the lightest Higgs boson around the 
LEP limit is considered to be more difficult. It necessitates more statistics 
(around 30 fb~^) than for supersymmetry, where 1 fb~^ are sufficient to be 
sensitive to masses of squarks of several hundred GeV/c^. The complexity 
in the Higgs sector is due to the necessity to search for rare, but clean fi- 
nal states (77) in large QCD/Standard Model backgrounds. Once the new 
particles are discovered, detailed measurements can be performed. 

To reconstruct the fundamental parameters, precise theoretical predic- 
tions are necessary. An overview of the tools is given in Ref . [HI] • The pro- 
duction rates for the Higgs boson are taken from fSU] . The most important in 
the supersymmetric sector are SuSpect, SOFTSUSY and SPheno [62 |[63l[M] 
for the calculation of the sparticle and Higgs masses. SUSY-Hit provides 
the branching ratios via SDecay for the decay of supersymmetric particles 
and HDecay for the decays of the Higgs bosons [651 EHl [67]. e"'"e~ cross 
sections are calculated by PYTHIA and SPheno [HSl El], while proton- 
proton cross sections, first calculated in |83j, are provided by Prospino2.0 
at NLO [MlllQlinilli. 

The pioneering work on extracting the fundamental parameters was per- 
formed in [73l [71], followed by Fittino [75] and SFitter [35]- The package 
Super-Bayes originated in the study of the dark matter aspect of supersym- 
metry [TB]. Several other projects are also underway. 

In this chapter, first the signatures of the Standard Model Higgs boson 
sector and the supersymmetric signals will be discussed. Then the most fa- 
vorable case will be studied: mSUGRA including all supersymmetric signals, 
followed by the study of the weak scale model, the MSSM. As a final study 
the question of the determination of the Higgs couplings as well as the sen- 
sitivity of these to new physics will be addressed. For the latter study the 
question will be: if new physics signals are not observed at the LHC, will the 
precision of the measurement of the Higgs boson couplings be sufficient to 
be sensitive to new physics? 

4.1 Low Mass Higgs Boson Signatures 

The Higgs boson can be produced in gluon fusion [85] , weak boson fusion and 
top pair production. The expected cross sections are shown in Figure 14.11 
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Figure 4.1: The production cross sec- 
tion for a Standard Model Higgs bo- 
son at the LHC is shown as function 
of its mass for the main production 
mechanisms. 
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Figure 4.2: The branching ratios of 
the Standard Model Higgs boson are 
shown as function of its mass. 



While tremendous progress on the understanding of the Monte Carlo mod- 
eling of the background processes has been made, it is essential for most of 
the analyzes to determine the background from the data, e.g., via side-band 
background subtraction. This means that at the LHC channels with the 
capability of the reconstruction of the Higgs mass (with the best possible 
precision) are favored over counting analyses. 

The gluon fusion process at the LHC is known with an uncertainty of 
13% [ng [1371 [1381 dSHl [HO]. NLO corrections [1301 [IM] as well as NNLO 
calculations with an effective ggH coupling |132[ I133[ I134[ 1135] lead to a 
cross section of 37 pb for a 120 GeV/c^ Higgs boson. For the weak boson 
fusion process, the cross section is about 4.5 fb |141l 11421 11431 1144] . This 
production mode relies on a central jet veto between the forward tagging 
jets |1451 11461 1147] . For this signal, NLO rates as well as distributions have 
been calculated [UHl EM EM [iMl [152] and an error of 7% [UHl EM EM 
[Tsn [T52] is used. 

The production channels of gluon fusion and weak boson fusion are not 
independent: Higgs boson production with two jets has contributions from 
both gluon fusion and weak boson fusion. With appropriate cuts, one channel 
can be enhanced over the other |1581 1159] . The error on this classification is 
covered by the error estimate. 

A cross section of 450 pb is expected for the channel ttH, the radiation 
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Table 4.1: Signatures used for the analysis for a Higgs mass of 120 GeV/c^. 
The Standard Model event numbers for 30 fb~ include cuts [86]. The factor 
after the background rates describes how many events are used to extrapolate 
into the signal region. The last two columns give the one-sigma experimental 
and theory error bars on the signal (taken from [36]). 

of a Higgs boson from a top quark. NLO rate calculations have been per- 
formed |160[ I16H 1162] . A theory error of around 13% is associated to this 
channel. 

For the associated production of the Higgs boson and a weak vector boson 
(Z or W), a cross section of 2.2 pb is expected. For this channel a theory 
error of 7% is assumed. 

The analysis will be constrained to the single production of Higgs bosons. 
Pair production is essentially sensitive to the self-coupling of the Higgs and 
extremely difficult at the LHC [1631 MM MB UMl HEZl MSI UMl EM [m] . 

For a Higgs boson of mass 120 GeV/c^ the highest branching ratio is 
to a pair of b-quarks. The branching ratio is roughly 90% as shown in 
Figure 14.21 However, discovering the Higgs boson in this decay channel is 
particularly difficult, because reconstructing a Higgs boson invariant mass 
peak from two b-jets above the QCD continuum is hopeless. Therefore the 
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associated production in combination with a pair of top quarks have been 
studied. The combinatorial background is more complex than previously 
estimated, making it difficult to establish this signal at the 5a level |171 [31] 
for 30 fb~ . This development is taken into account by reducing the signal 
by 50% with respect to the analysis in Ref. [86] . 
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Table 4.2: Experimental systematic errors used in the analysis. Left: sys- 
tematic errors applying to both signal and background. Reconstruction and 
tagging efficiencies are defined per particle, e.g., ii —)■ 77 has a 4% error on 
the photon reconstruction. Right: systematic background errors, either fully 
correlated or independent between channels. Tables are the same as [86tl36]. 

The decay of the Higgs boson to two photons has the advantage of a pre- 
cise mass reconstruction. Thus in spite of the low branching ratio of order 
of 10~^, the excellent mass resolution of the calorimeters at the LHC allows 
to separate the signal from the background of QCD and prompt photons. 
As the natural width of the Higgs boson at 120 GeV/c^ is smaller than the 
experimental resolution (about 2 GeV/c^), the uniformity of the calorime- 
ters is essential for the discovery. For the measurement of the Higgs boson 
mass, once sufficient statistics are accumulated, a precise knowledge of the 
photon energy scale in the calorimeter is essential. The preparation of this 
measurement through test beam studies gives additional confidence in this 
channel. The channel can be measured separately in the gluon fusion, weak 
boson fusion and top pair production. The expected precision of the mass 
measurement is (9(200MeV/c^). 

The associated production of ZH or WH was considered to be hopeless 
initially (in contrast to the TeVatron where it is the main source for the 
signal search). In a theoretical study the Higgs boson is required to be 
strongly boosted with an inherent loss of efficiency due to this requirement. 
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The idea is then to use the subjet structure of the Higgs boson decay to gain 
additional leverage |SZ] more than compensating the loss of efficiency. This 
promising new channel is also included in the analysis. The results of the 
theoretical analysis have been confirmed recently by an ATLAS study at the 
level of 10% Hg. 

The branching ratio for Higgs boson to a pair of Z bosons is smaller at 
120 GeV/c^ than at much higher masses. It is considered to be the golden 
channel with the Z bosons decaying to electrons or muons. Some information 
can still be obtained from this channel, either produced in gluon fusion or in 
weak boson fusion. 

The decay of the Higgs boson to two leptonically (electrons and muons) 
decaying W^ bosons is also used in the analysis. This counting analysis, due 
to the neutrinos in the W^ decay, does not allow the reconstruction of the 
Higgs boson mass. The background is rejected using the angular correlation 
of the leptons. 

The decay of the Higgs boson to a pair of r-leptons can be used for the 
weak boson fusion process. Here the Higgs boson is boosted (in contrast to 
the gluon fusion process) which allows the use of the collinear approximation 
to reconstruct the Higgs boson mass. 
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Table 4.3: Theory errors used in the analysis. The left table shows the errors 
on the production cross section and the right table shows the errors on the 
branching ratio (taken from [36]). 

The production rates for the Higgs sector were calculated in [HH] and the 
branching ratios are taken from HDecay |67[[66]. All rates and experimental 
errors are effectively taken from [86] (with the exception of the modified 
ttbb channel). The expected results are listed in Table WA\ Experimental 
systematic errors are shown in Table W2\ and the theoretical errors are shown 
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in Table E3\ 

At the ILC the mass of of a 120 GeV/c^ Higgs boson can be measured 
with a precision of about 50 MeV/c^. The expected precision on the Higgs 
couphng measurements is at the percent level. Additionally precision mea- 
surements of the couplings of the Higgs boson are not restricted to judiciously 
chosen channels. A model independent analysis of the Higgs production and 
decay is possible. This can have consequences on the interpretation of the 
determination of the couplings as will be shown later. 

4.2 Super symmetric Observables for a spe- 
cific Parameter Set 
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Figure 4.3: The production cross sec- 
tion for supersymmetric particles at 
the LHC at NLO is shown as calcu- 
lated by Prospino2.0. 
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Figure 4.4: The effective mass is 
shown for the Standard Model and 
several supersymmetric parameter 
sets in ATLAS. 



The parameter point SPSla is characterized by moderately heavy squarks 
and gluinos, which leads to long cascades including neutralinos and sleptons. 
The mass of the lightest Higgs boson is close to the mass limit determined 
at LEP. 

The production of squarks and gluinos and their decay lead to large sig- 
nals (in contrast to the Higgs sector) above the Standard Model background. 
The NLO cross sections calculated by Prospino2.0 are shown in Figure 14.31 
The production cross sections for strongly interacting particles are several 
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tens of picobarn. Inclusive analyses such as the search for highly spherical 
events with large missing transverse energy (due to the LSPs) and several 
jets with high transverse momenta, and/or adding the requirement of a lep- 
ton (or more) will allow the discovery of supersymmetry. The effective mass, 
the sum of the jet transverse momenta and the missing transverse energy, 
will give a first hint of the scale of supersymmetry. An example for the ef- 
fective mass of supersymmetric signals and Standard Model background is 
shown in Figure 14.41 At large effective mass the supersymmetric signal dom- 
inates the distribution. It is essential to have a good understanding of the 
jet and ETmiss reconstruction to have confidence in the discovery. In the 
standard scenarios the LHC will be sensitive to squark and gluino masses up 
to 2.5 TeV/c2. 

An important aspect of subsequent measurements of the supersymmetric 
particles is the long decay chain |106[ 1107] 

qL ^ xk ^ ^B.^(i ^ (^(^qxv (4.1) 

The final state contains at least a hard jet and two opposite sign same flavor 
leptons. In this decay chain five edges and thresholds can be calculated and 
reconstructed |16]. One of the invariant mass combinations is the lepton- 
lepton edge. As shown in Table H^ with an integrated luminosity of 300 fb~^, 
the measurement is not statistics limited, but limited by the knowledge of 
the lepton energy scale. Thus attaining the goal of a per mil level energy 
scale is important and justifies the tremendous work performed on this issue 
in test beam and Monte Carlo simulations. 

The edges and thresholds can be expressed as functions of the four in- 
tervening masses. Therefore from the system shown in Table 14.41 the masses 
of the supersymmetric particles can be reconstructed, using either via toy 
Monte Carlo or a fit. It is important to note that there is no assumption is 
made on the underlying theory. Further signatures, e.g. the squark-R and 
the sbottoms, provide a total of 14 observable and measurable particles at 
the LHC in SPSla. The precision of the mass determination for the LHC 
is shown in Table 14.51 Typically the systematic error on measurements at 
the LHC coming from the jet energy scale is 1% and 0.1% for the lepton 
energy scale. These energy-scale errors are each taken to be 99% correlated 
separately as discussed in Ref [46j. With integrated luminosities of up to 
300 fb~^ the statistical error in many cases is smaller than the systematic 
error. 
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Table 4.4: LHC measurements in SPSla, taken from |16]. Shown are the 
nominal values (from SuSpect) and statistical errors, systematic errors from 
the lepton (LES) and jet energy scale (JES) and theoretical errors. All values 
are given in GeV. Table taken from [35] . 



For the ILC, as a rule of thumb, if particles are light enough to be pro- 
duced in pairs given the center-of-mass energy of the collider, their mass can 
be determined with an order of magnitude better precision. Two separate 
methods provide mass measurements. Direct reconstruction, i.e., for the ob- 
served decay products of the sparticles in the detector relies on the precision 
of the energy scale of the detector. The second method, a counting exper- 
iment, is based on the measurement of the cross section as function of the 
(variable) center-of-mass energy of the ILC. Here, experimentally, the pre- 
cise determination of the beam energy is of the utmost importance, whereas 
the challenge for the detector data is the understanding of the efficiency. 
Both methods have comparable precision in the per mil region, but different 
systematical errors. Additionally discrete quantum numbers like the spin 
of the particles can be determined as well by the analysis of the angular 
distributions. 
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Table 4.5: Errors for the mass determination in SPSla, taken from ^6] . 
Shown are the nominal parameter values (from SuSpect), the error for the 
LHC alone, from the ILC alone, and from a combined LHC+ILC analysis 
(taken from [35] ) 



The summary of particle mass measurements is listed in Table 14.51 taken 
from Ref. |35]. The central values are calculated by SuSpect |62]. The 
coverage of the LHC is better in the strongly interacting sector, whereas the 
ILC has the advantage of higher precision in the gaugino/slepton sector [901 
Ion [921 [93] . As pointed out in Ref. [721 [3]) the use of the ILC measurement 
of the LSP in the kinematic determination of the squark masses increases 
their precision. 



4.3 Treatment of Errors 

A rigorous treatment of the errors is essential to obtain the most precise 
and also reliable estimation of the fundamental parameters. There are three 
different sources of errors: experimental statistical errors, experimental sys- 
tematic errors and theoretical errors. Systematic errors can have different 
sources, such as the energy scale or the knowledge of the efficiency and vary- 
ing degrees of correlation. The picture even depends on the measurements 
used in the analysis: edges and thresholds have independent statistical errors 
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and correlated energy scale errors, but when translating them into masses, 
the masses are then non-trivially correlated. Theoretical errors are an esti- 
mate of unknown higher order corrections. 

In SFitter the CKMfitter prescription for the RFit scheme [HU |95] is 
followed. The experimental errors are treated as Gaussian and the full cor- 
relation matrix is included. The theoretical errors are fiat errors. Therefore 
any central value within the range is equally probable. This description is 
particularly well adapted as the error covers also missing higher order cal- 
culations. In this case the higher order, once it is calculated, should move 
the central value of the prediction within the error. This behavior would be 
contrary to the Gaussian expectation. 

While in the supersymmetric measurements the errors are essentially 
Gaussian, for the study of the Higgs parameters, low counting rates (af- 
ter background subtraction) are expected. A combination scheme of Poisson 
and Gaussian errors was developed under the guiding principle that in the 
Gaussian limit the formula should simplify to the quadratic sum of the errors. 
The detailed procedure is described in the Appendix of [5B] . 

The contribution to the x^ of a given measurement is zero within la of 
the theoretical error and outside of this range the experimental error is used. 
Given a set of measurements d and a general correlation matrix C 



X^ = Xd C ^ Xd 



\dj — dA < (7, 



(theo) 



\Xd,i\ = { M.-d-i-^f"""' I , _ J I . ^(theo) (4.2) 

i 

where di is the i-th data point predicted by the model parameters and di 
the measurement. 

A complication with fiat distributions is that in the central region the 
log-likelihood can be constant as a function of model parameters. In those 
regions these parameters vanish from the counting of degrees of freedom. 

To determine the errors on the fundamental parameters, two techniques 
are used: a direct determination for the best fit using MINUIT and using 
sets of toy measurements. The advantage of MINUIT is that only one fit is 
necessary to determine the errors, but the presence of correlations and the 
fiat theory errors complicate the error determination. The toy experiments 
method is used by SFitter unless stated otherwise. The search for the best-fit 
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point is performed several (10000) times. The datasets are smeared according 
to the experimental and theoretical errors including the correlations among 
them. 

4.4 Measuring mSUGRA 

Determining the parameters from a given set of measurements necessitates 
efficient algorithms. The determination of the mSUGRA parameters can be 
used as an illustrative example. Four parameters have to be determined and 
one sign. The smuon mass depends on uiq and mi/2, the neutralino mass 
depends on tan/3 and mi/2, among others. The measurements are strongly 
correlated, e.g., via the energy scale error. 

4.4.1 Determination of the True Central Value 

The first question to be asked is whether the correct parameter set can be 
found. Thus the algorithm has to be able to disentangle the different con- 
tributions. The simplest approach is to use a gradient fit, i.e., MINUIT. In 
a tightly constrained system such as mSUGRA this is sufficient to find the 
correct parameter set even if the starting values for the parameters are far 
from the true SPSla values. The sign of the n parameter is then determined 
by the quality of the fit: performing the fit with the wrong sign of /x leads to 
larger x^ value than the correct sign of /i. 

However, the result could depend on the starting point when the fit re- 
mains confined to a secondary minima by a "potential" wall in x^- Simulated 
annealing was developed in Ref. ^5j to allow to escape the secondary mini- 
mum and find the correct solution. In SFitter weighted Markov chains were 
developed. Markov chains are linear in the number of parameters, an im- 
portant property for the determination of the MSSM parameters, and are 
capable of finding secondary minima. 

First, fully exclusive log-likelihood maps are calculated. From these dif- 
ferent types of projections are possible: Bayesian or frequentist (profile like- 
lihood). Successive projections over unwanted parameters allow to study 
the correlations between parameters. In the Bayesian approach a measure, 
the Bayesian prior, is introduced and the parameter is integrated over. In 
the frequentist approach the parameter value of the minimum x^ (maximum 
likelihood) is retained. The Bayesian approach ensures that after the pro- 
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jection the distributions have the quahty of a probabihty density function. 
The advantage of the frequentist approach is that the absolute minimum is 
always retained, however the pdf property is lost. 
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Figure 4.5: mo-mi/2 plane of SPSla Figure 4.6: mo-mi/2 plane of SPSla 
for LHC measurements after the for LHC measurements after the fre- 
Bayesian projection (from [35]). quentist projection (from [35J). 



An example is shown in Figure 14.51 for the Bayesian approach with a flat 
prior and in Figure 14.61 for the frequentist approach in the mo-mi/2 plane 
of SPSla for LHC measurements. While the general structure is similar 
and the correct parameters were determined, due to the definition of the 
frequentist projection the best-fit is brighter. In the Bayesian case the sharp 
maximum/minimum is washed out by the integration over large regions with 
small likelihood (noise effect). One can also see the parabolic correlation 
between the parameters as expected from, e.g., the sfermion masses. 
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Table 4.6: List of minima found in mSUGRA for the LHC measurements 
(taken from [33]). 

The expected error on the top quark mass measurement is 1 GeV/c^ at the 
LHC and an order of magnitude more precise at the ILC (0.12 GeV/c^). The 
top quark mass therefore has a double role, it is a parameter of the model [96| 
|97] [98] IM] which can have a large impact on the prediction of masses such 
as the Higgs boson mass [HSl llOOj 1101] . The top quark mass is also treated 
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as a measurement with the appropriate error in the determination of the 
parameters. In principal all Standard Model parameters would need to be 
treated this way, however, the impact is practically negligible for the other 
parameters, so that they can be fixed to the central value. An exception is 
the b-quark mass for high tan /3 regions. 

The list of secondary minima for the LHC measurements is shown in 
Table 1121 including the top quark mass as a parameter. Even in such a tightly 
constrained model with a small number of parameters secondary minima can 
be observed in addition to the correct parameter point (first line of Table S]6]). 
It is particularly interesting to note that the addition of the top quark mass as 
a parameter adds a secondary minimum (third line) in its interplay with tri- 
linear coupling. However, all secondary minima can easily eliminated by the 
absolute value of the x^ which is higher than for the correct SPSla parameter 
set. A dataset which was smeared according to the expected errors, including 
all correlations, was used for this determination. Therefore a small shift in 
the parameters is expected as the dataset does not correspond to the true 
value of SPSla. 

4.4.2 Determining the Errors on the Parameters 

The second question is to determine the precision with which the parameters 
can be determined. At this point the theoretical errors have to be specified. 
The Higgs boson mass prediction is known with a precision of 3GeV/c^. 
Strongly interacting particles such as squarks and gluinos have an error of 3%. 
Electromagnetically and weakly interacting particles such as neutralinos and 
sleptons have an uncertainty of 1% of their mass prediction [9T| [921 193]. The 
errors defined on the masses were propagated to the edges and thresholds. 

At the LHC, one has the choice to either use the edges and thresholds 
directly or the masses derived from the edges and thresholds |106[ I107[ I108[ 
1109] . A comparison of the expected precision is shown in Table 14.71 The pre- 
cision is improved by the use of edges and thresholds instead of the masses. 
In fact, the determination of the mass introduces non-negligible correlations 
among them. As these correlations are not available, there is a loss of infor- 
mation which translates into a less precise determination of the fundamental 
parameters by a factor of between two and three for mo and mi/2. Further 
evidence for the impact of correlations is seen by the improvement of the 
error on uiq when the correlation of the lepton and jet energy scale errors are 
taken into account. 
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Table 4.7: Errors on the parameter determination at the LHC in SPSla. 
The big columns are mass and endpoint measurements. The subscript repre- 
sents neglected, (probably approximate) Gaussian or flat theory errors. The 
experimental error includes correlations unless indicated otherwise in the su- 
perscript (taken from [35]). 

Theory errors have an impact on the determination of the parameters at 
the LHC. As shown in Table 14. 7[ the proper inclusion of theory errors as 
fiat errors reduces the LHC precision significantly. An example is the Higgs 
boson mass which is sensitive to tan/3, but its theoretical error is an order 
of magnitude worse than the expected experimental error. 

Theoretical errors not only parametrize not-yet calculated high order cor- 
rections but also variations of results between calculations at the same order. 
As a checlc of the consistency of the theory errors, the SPSla spectrum cal- 
culated by SoftSUSY was fitted with SuSpect. While the parameter central 
values were shifted, they were within 3a of the correct ones adding confidence 
in the theory errors. 

The table of measurements at the LHC hides a part of the complexity 
of the determination of the parameters: three neutralinos are observed, the 
third lightest, due to Higgsino nature is not. Therefore one has to be able to 
infer from the data which of the neutralinos were observed. The most difficult 
example is to wrongly label the measurement of the fourth neutralino as third. 
A parameter determination can be performed and mg and mi/2 are off by up 
to 1 GeV/c^. The x^ value of the fit remains reasonable. While this looks 
like a true difficult secondary minimum, the parameters determined can be 
used to predict the spectrum, production cross section and branching ratios. 
Analyzing these results, a larger production of Xa than Xs is expected in 
disagreement to the labeling, discarding effectively this secondary minimum. 
It is clear that a lot of thought has to go into the definition of the table of 
observables to make sure that all underlying/hidden assumptions are tested. 
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An automated algorithm for all does not seem feasible. 
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Table 4.8: Expected errors for mSUGRA at the LHC (endpoints) and ILC 
Only absolute errors are given. Flat theory errors are used (taken from [53]). 



The ILC will start taking data after the LHC has taken data and dis- 
covered SPSla (if nature is kind enough). The added value of the ILC is 
the completion of the particle spectrum and the measurement of all particles 
with an impressive accuracy, about an order of magnitude improved on the 
LHC. The results are shown in Table WM The precision improvement by the 
ILC per se is about a factor of five. The improvement of the combination 
LHC and ILC over ILC alone is rather marginal. Comparing the LHC+ILC 
errors with and without theory errors shows the margin for the improvement 
of theory predictions, as shown already for the LHC, justifying the SPA 
project |13]. The SPA project's goal is to improve the theory predictions for 
the supersymmetry (S) parameter (P) analysis (A). 

4.5 Measuring the MSSM 

mSUGRA has the advantage of being a model with few parameters which 
is mostly defined at the GUT scale. The unification of breaking parame- 
ters at the GUT scale is assumed. The MSSM, defined at the electroweak 
scale, allows to measure grand unification without imposing it. The MSSM 
parameter determination is far more difficult technically: more parameters 
have to be determined from the same dataset. The reward is to be able to 
extrapolate the MSSM parameters to the GUT scale and thus measure grand 
unification or other models such as GMSB. 
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Table 4.9: Eight best-fitting points in the MSSM with two ahernative solu- 
tions for At. The x^ value for all points is approximately the same, so the 
ordering of the table is arbitrary. The parameter point closest to the correct 
point is labeled as SPSla (taken from [35j). 
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4.5.1 Determination of the True Central Value 

The problem of the higher-dimensional MSSM parameter space with respect 
to the mSUGRA case is solved by a combination of techniques (for details 
see [35]). In particular a sequence of Markov chains and gradient fits is used. 
A second Markov chain step combined with the gradient fit is performed in 
the gaugino-Higgsino subspace and the other parameters are re-determined 
subsequently. This complex procedure leads to the identification of the pri- 
mary minimum/maximum and an eightfold ambiguity listed in Table 14.91 
Additionally, alternative likelihood maxima are triggered by correlations be- 
tween the rather poorly measured parameters At, tan/3 and the right-handed 
stop mass. 
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Figure 4.7: The Bayesian projection Figure 4.8: The profile likelihood for 

for Ml is shown for LHC data [35]). Mi is shown for LHC data [35]). 

Solutions are observed for the two signs of fi and for each sign of /i four 
solutions are determined, two of which can be grouped together by inter- 
changing Ml and M2. The correct solution differs from the nominal values 
of SPSla as a smeared dataset was used for this study. In Figure 14.71 the 
Bayesian projection is shown and in Figure H^ the frequentist projection is 
shown. While in the Bayesian projection the central values are washed out 
by noise effects, the profile likelihood allows to identify the location of the 
ambiguities. 

4.5.2 Determining the Errors 

In principle 22 measurements are available at the LHC, thus technically 
enough to determine the 19 parameters of the MSSM (including the top 
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quark mass as a parameter). However, these measurements come from only 
15 different masses. Therefore, inevitably, some parameters must be fixed in 
order to determine the errors with no impact on the determination as these 
parameters have essentially a very weak impact on the observables. The 
result is shown in Table 14.101 including fiat theory errors. 

The stau sector is inevitably poorly determined at the LHC as only one 
measurement is available for the parameters. The stop sector (no measure- 
ments in SPSla) can only be determined indirectly with a penalty paid in 
the size of the error. For the ILC several parameters have to be fixed as the 
strongly interacting particles such as squarks and gluinos with the exception 
of the lightest stop cannot be produced. The historical motivation for fixing 
the squark parameters in the ILC study was the possibility of concurrent 
running of LHC and ILC. As this is no longer an option the more realistic 
scenario will be that the parameters parameters for the ILC-alone study will 
be fixed to the measurements of the LHC. The fixed parameters in the ILC 
column of Table 14.101 are within 5%-15% of such a scenario. The main im- 
pact of the combination of LHC and ILC is observed in the last column of 
Table 14.10^ where all parameters can be determined, with the exception of 
the tri-linear couplings of the first two generation, which have no effect on 
the observables in the selectron and smuon sector. 

Comparing the effect of taking into account the theoretical errors, about 
a factor five in precision is lost with today's calculations at the ILC. Its 
precision is decreased from per mil to 0.5%. For the LHC roughly a factor two 
is lost. The combination of LHC and ILC measurements can be particularly 
useful to determine the link to dark-matter observables [791 1122[ I123[ I124[ 
[T25| [76| fT26| [127 1 [128]. 

4.5.3 Extrapolating to the GUT Scale 

Once the parameters of the weak-scale MSSM-Lagrangian have been de- 
termined, the next step is to extrapolate the parameters all the way to the 
Planck scale. Inspired by the apparent unification of the gauge couplings [T7] 
in the MSSM the question arises if any other running parameters unify at a 
higher scale, as shown in the pioneering work in [731 El] • 

Technically, upwards running is considerably more complicated |129] than 
starting from a unification-scale and testing the unification hypothesis by 
comparing to the weak-scale particle spectrum. For example, it is not guar- 
anteed that the renormalization group running will converge for weak-scale 
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199.6±C(104) 


-763.4 


All,2 


fixed 


fixed 


fixed 


-251.1 


Aul,2 


fixed 


fixed 


fixed 


-657.2 


Adl,2 


fixed 


fixed 


fixed 


-821.8 


rriA 


406.3±C(103) 


393.8± 1.6 


393. 7± 1.6 


394.9 


fi 


350.5± 14.5 


354.8± 3.1 


354.7± 3.0 


353.7 


ITLt 


171.4± 1.0 


171.4± 0.12 


171.4± 0.12 


171.4 



Table 4.10: Results for the general MSSM parameter determination in SPSla 
using fiat theory errors. The kinematic endpoint measurements are used for 
the LHC and the mass measurements for the ILC. The LHC-I-ILC column 
is the combination of the two measurements sets. Shown are the nominal 
parameter values and the result after fits to the different data sets (taken 
from [351. 
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Figure 4.9: Extrapolation of the in- 
verse of the gaugino mass parameters 
to the GUT scale for the correct so- 
lution at the LHC 




Figure 4.10: Extrapolation of the in- 
verse of the gaugino mass parameters 
to the GUT scale for a false solution 
at the LHC. 



input values far away from the top-down prediction. 
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Figure 4.11: Extrapolation of the 
first and second generation scalar 
mass parameters to the GUT scale 
for the correct solution at the LHC. 



Figure 4.12: Extrapolation of the 
first and second generation scalar 
mass parameters to the GUT scale 
for a false solution at the LHC. 



The first step in the extrapolation of the weak scale parameters is to 
test the expected measurement of the gaugino mass parameters at the LHC. 
Figure 113] shows the result for the true parameter point. In Figure H. 101 the 
extrapolation is performed for one of the ambiguous solutions. While grand 
unification is observed for the correct parameter set as expected, in the second 
case unification is not observed. The non-unification of six of the solutions is 
easy to observe. The seventh solution shows a close miss of unification. It will 
therefore depend on the exact precision of the LHC measurement whether it 
can be distinguished from the correct solution which unifies the parameters 
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at the GUT scale. 

For the extrapolation of the scalar mass parameters at the LHC alone the 
picture is even more difficult. In Figure H. Ill the unification of the first and 
second generation scalar mass parameters is observed for the correct solution, 
however in Figure 14.121 the stau parameters have been moved to different 
values. This behavior can be understood from the RGEs. In the gaugino 
sector the RGEs do not depend on other parameters. In the scalar sector the 
RGEs are strongly coupled, so that a wrong value of a stau parameter can 
wreak havoc in the extrapolation of the first and second generation scalar 
parameters. 

If no further measurements are thought of for the SPSla scenario, the 
proof and measurement of grand unification of supersymmetric parameters 
will have to wait until the ILC resolves the ambiguities. With the ILC a 
unique solution of the MSSM parameters is determined, which can then 
be extrapolated to truly measure grand unification of the supersymmetric 
breaking parameters as shown in Chapter |5l 

4.6 Measuring the Higgs Sector 

A more difficult scenario at the LHC would be to discover only a Higgs 
boson in two-particle decays. The question arises naturally if there is more 
to learn from this sector than the Higgs boson mass, which is used also in the 
supersymmetric parameter determination. In analogy to the measurement of 
the Triple Gauge Couplings (TGC) at LEP, one can parametrize the Higgs 
couplings as deviation from the Standard Model expectation. 

Thus a tree-level Standard Model Higgs coupling to any particle j is 
defined as: 

am-^g^j^H a + ^m) (4-3) 

where the Ajjh are independent of each other. The three loop induced cou- 
plings in the Higgs sector, QggH, g-yyu and g^zH, sensitive to new physics |172[ 
I173[ I174[ 1175] , are defined in the following way: 

gm -^ g^H (1 + ^!h + A,,^) (4.4) 

where g^^ is the loop-induced coupling in the Standard Model, A*^*^^ is the 
contribution from modified tree-level couplings to Standard-Model particles 
(in the loop), ^jjh is an additional contribution proper for this effective 
coupling. 
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In the scenario treated in this analysis, a Higgs boson mass of 120 GeV/c^, 
the total width is defined as the sum of the measured partial widths. As the 
partial width of the decay Higgs to bb is 90% of the total width, it is measured 
reasonably well at the LHC 

The total width can be varied just like the couplings 

Ttot = rs„^ (1 + Ar) . (4.5) 

This parameter can be used of an unobserved channel/anomalous channel. 

4.6.1 Determination of the True Couplings 




HWW 



^ -1 i 



Figure 4.13: Profile likelihood for the 
ttH and WWH couplings for 30 fb"^ 
(from [36]). 



^ 




-10 12 3 

HWW 

Figure 4.14: Bayesian projection for 
the ttH and WWH couplings for 
30 fb~^ (from p]). 



The procedure for the determination of the Higgs boson couplings fol- 
lows the strategy applied in the determination of the supersymmetric cou- 
plings. An exclusive likelihood map is calculated for the parameters AjjH, 
the deviations from the Standard Model couplings. From this map different 
projections are possible to study correlations. 

The analysis is performed for an integrated luminosity of 30 fb^ . The 
errors cannot be scaled simply to 300 fb~ since e.g., the systematic error 
on the mini-jet veto for nominal LHC luminosity is not known. The first 
parameter set used for the true dataset (without smearing) are the tree-level 
couplings: WWH, ZZH, ttH, ttH and no genuine anomalous contribution 
to the effective couplings is allowed. The overall phase of the couplings is 
fixed by requiring WWH to be positive without loss of generality. 
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While in the supersymmetric case genuine secondary minima can arise, 
e.g., particularly in the MSSM, in the Higgs couplings no such effect is ex- 
pected. In fact, the definition of the parameters leads to quadratic depen- 
dence of the observables on the parameters. 

The Bayesian projection is shown in Figure 14.131 and the frequentist ap- 
proach of maintaining always the absolute minimum/maximum is shown in 
Figure 14.141 From the comparison of the two figures one can conclude that 
the important integration (noise) effect renders the Bayesian approach less 
performant than the profile likelihood. 

The profile likelihood in Figure 14.141 shows a well identified correct solu- 
tion. It is interesting to note that in spite of the observables being essentially 
quadratic in the parameters, a clear preference for the correct solution is ob- 
served for Htt. The reason is that the rate of Higgs to two photons is sensitive 
to the relative sign. 

In general, a positive correlation among the couplings is observed among 
the couplings induced by the general (a ■ BR) structure of the LHC measure- 
ments. A change in the bbH coupling will induce a coherent compensation 
on the production side and decay side of the observables as the total width 
in the denominator is dominated by the bbH partial width. 







HWW 




Figure 4.15: Profile likelihoods for 
WWH versus ttH allowing genuine 
contributions to the effective cou- 
plings (taken from [36]). 
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HWW 

Figure 4.16: Profile likelihoods for 
WWH versus bbH allowing genuine 
contributions to the effective cou- 
plings (taken from 



In Figure 14.151 and Figure 14.161 the correlations of the WWH with the 
ttH and bbH couplings are shown. Here the effective couplings are left free 
and are part of the parameter set to be determined. The general feature 
of a positive correlation is still observed, especially clear for the bbH versus 
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WWH coupling. However, the sign sensitivity of the 77 channel is lost due 
to the additional freedom added by the effective anomalous couplings. 

While these results seem to indicate that the absolute couplings of the 
Higgs boson can be determined at the LHC, a simple scenario shows that this 
is not the case. In the studies discussed above, the ccH coupling was fixed 
to its Standard Model value. What would be the effect of a large anomalous 
ccH coupling? The observed rates would be diminished by Th/Th+ccH- In 
the coupling determination the anomalous couplings would deviate from the 
Standard Model value and be negative. This situation is thus indistinguish- 
able from genuine anomalous couplings. 

In practice the coherent observation of the reduction of all couplings would 
lead to a study of the Ap parameter. This parameter is determined in a single 
parameter fit. It cannot be determined in a global fit including the anomalous 
Higgs couplings due to the structure of the observables (blind direction). 

Absolute coupling measurements at the LHC therefore depend on the 
(additional) assumptions to be made in contrast to coupling ratios. A model 
independent study will have to wait for the ILC 

4.6.2 Error determination 





no effective coupling 


s 


with effective couphn 


gs 




RMS 


'^symm 


O^neg 


'^pos 


RMS 


O^symm 


'^ncg 


'^pos 


TUh 


±0.36 


±0.26 


-0.26 


±0.26 


±0.38 


±0.25 


-0.26 


±0.25 


^WWH 


±0.31 


±0.23 


-0.21 


±0.26 


±0.29 


±0.24 


-0.21 


±0.27 


^ZZH 


±0.49 


±0.36 


-0.40 


±0.35 


±0.46 


±0.31 


-0.35 


+ 0.29 


AttH 


±0.58 


±0.41 


-0.37 


±0.45 


±0.59 


±0.53 


-0.65 


+ 0.43 


\-bH 


±0.53 


±0.45 


-0.33 


±0.56 


±0.64 


±0.44 


-0.30 


+ 0.59 


ArfH 


±0.47 


±0.33 


-0.21 


±0.46 


±0.57 


±0.31 


-0.19 


+ 0.46 


A-y-yH 


— 


— 


— 


— 


±0.55 


±0.31 


-0.30 


+ 0.33 


A99H 


— 


— 


— 


— 


±0.80 


±0.61 


-0.59 


+ 0.62 


nib 


±0.073 


±0.071 


-0.071 


±0.071 


±0.070 


±0.071 


-0.071 


+ 0.072 


rrit 


±1.99 


±1.00 


-1.03 


±0.98 


±1.99 


±0.99 


-1.00 


+ 0.98 



Table 4.11: Errors on the measurements from 10000 toy experiments for an 
integrated luminosity of 30 fb~^ (taken from [36j ) . 

The errors on the Higgs boson couplings are determined with 10000 toy 
datasets. The datasets were smeared according to the experimental and 
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theoretical errors, including the correlations among them. The results are 
shown in Table 14.111 without and with anomalous effective couplings. 

The RMS is shown, as well as the fit of a symmetric Gaussian and the fit of 
two separate Gaussian for the positive and negative side of the coupling. The 
expected precision for the coupling measurement is typically between 25% 
and 45%. As statistics increase (higher luminosity) the difference between 
RMS and the Gaussian fits should decrease. 

Adding the effective couplings to the parameter sets, one would expect an 
increase of the errors. While this is true in most cases, e.g., the ZZH coupling 
is decreased. One can understand this effect by the impact of correlations. 
The rate of the ZZ final state is the upper limit. With a positive correlation 
between the effective coupling (production side) and the HZZ (decay side), 
the effective values that the ZZH coupling can take on are actually reduced, 
leading to a smaller error. 





no effective couplings 


with effective couplings 




'-'"symm 


'^neg '^pos 


'^symm 


(Tjieg 


'^pos 


^ZZhI ^WWH 


±0.46 


-0.36 ±0.53 


±0.41 


-0.40 


±0.41 


^tm/^wwH 


±0.30 


-0.27 ±0.32 


±0.51 


-0.54 


±0.48 


^hhu/ ^WWH 


±0.28 


-0.24 ±0.32 


±0.31 


-0.24 


±0.38 


^tth/ ^WWH 


±0.25 


-0.18 ±0.33 


±0.28 


-0.16 


±0.40 


/S.^^h/^wwh 


— 


— — 


±0.30 


-0.27 


±0.33 


'^ggn/^WWH 


— 


— — 


±0.61 


-0.71 


±0.46 



Table 4.12: Errors on the ratio of couplings, corresponding to Table 14.111 
(taken from [36j). 



The errors on the ratios are shown in Table 14.121 As many experimental 
and theoretical errors are correlated, one would expect that the ratios are 
measured more precisely than the absolute couplings. However the effect 
is relatively small. This is essentially due to the statistics of the signal for 
30 fb~^. For the experimental error the statistical error is much larger than 
the systematic error which is canceled in the ratio. The definition of the flat 
errors implies that outside the theoretical errors, the x^ varies only with the 
experimental error. As the experimental error is much larger, the theoretical 
errors do not have a large impact. This has been checked by removing the 
theoretical errors entirely, which has only a marginal effect. 

To illustrate the importance of the measurement of the bbH coupling, in 
Figure d. 171 the effect of the subjet analysis is shown with full sensitivity, re- 
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Figure 4.17: Impact of the subjet analysis. The sensitivity is 100% (left), 
50% (middle) and entirely removed (right) (taken from [36j). 

ducing its sensitivity by 50% and removing it entirely. ATLAS has confirmed 
the theoretical study within recently 10%, i.e., 3.7cr significance versus 4.2a 
which is closest to the full sensitivity in Figure 14.171 

4.6.3 Supersymmetric Higgs Scenarios 

In the favorable case of low mass supersymmetry, supersymmetric particles 
will be discovered well before the lightest Higgs boson. However, there are 
scenarios where the supersymmetric particles are out of reach or need very 
high integrated luminosity to be discovered and the Standard Model Higgs 
boson is the first "supersymmetric" particle discovered. In such a case it 
is interesting to ask whether the measurement of the Higgs boson couplings 
could be used to probe new physics in the absence of/ or ignoring direct 
evidence. 

Two scenarios were studied in [36j of which only one will be discussed 
here as an illustration. The parameter set is a modified version of SPSla. 
The study is based on the assumption that only a Higgs boson has been 
measured with a mass of 120 GeV/c^ and no supersymmetric measurements 
have been made. 





^WWH 


^ZZH 


Attn 


^bbH 


^rrH 


^-nH 


^ggH 


niH 


true 


-0.13 


-0.13 


-0.19 


3.27 


3.29 


0.19 


-0.28 


120.0 


errors 


±0.45 


±0.61 


±0.63 


±2.34 


±3.35 


±0.99 


±1.12 


±0.29 




-0.43 


-0.99 


-0.60 


-3.68 


-3.23 


-0.70 


-0.69 


-0.29 




+0.48 


±0.52 


±0.65 


±1.52 


±3.58 


±1.30 


±1.46 


±0.30 



Table 4.13: True couplings for the SPSla-inspired scenario. The errors were 
calculated for 30 fb~^ with 10000 toy experiments (taken from [36]). 
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The datasets were modified to take into account the modification of 
the couphngs. With respect to SPSla, tan/3 was modified to 7, At to 
— 1100 GeV/c^ and rriA to 151 GeV/c^. The true values for all Higgs cou- 
plings are listed in the first line of Table 14.131 Comparing these deviations 
with the Standard Model error bands in Table I4.1H it seems daunting to 
draw conclusions from the individual measurements, but in addition to the 
individual precision of the coupling measurements, the full correlation matrix 
(added information) is used. 

The first question to be asked is whether the Standard Model is excluded. 
The exclusion limit for the Standard Model is calculated by determining the 
X^ distribution for smeared Standard Model datasets, calculated with the 
fixed Standard Model couplings. The integral over x^ to 90% of the datasets 
gives the 90% confidence level. Then the new physics (SPSla-modified) 
datasets are used and the x^ is calculated with the Standard Model cou- 
plings. The two dimensional results are that for 77% of the toy experiments 
the Standard Model is excluded at a confidence level of 90%. A discovery 
is much more difficult, for 4% of the toy experiments the new physics sce- 
nario is a better description of the datasets than the Standard Model. Using 
the Higgs coupling measurements, one may be able to exclude the Standard 
Model, but a discovery will be difficult. 
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Chapter 5 
Conclusions 



After many years of preparation, the LHC has started operations in 2008. 
The first colhsions are expected in the run 2009/2010 which will start in 
November 2009. The increase of the center-of-mass energy beyond the TeVa- 
tron will open up a new energy regime to search for new physics. 

The first year(s) of LHC and ATLAS/CMS operations will concentrate 
on the understanding of the detector using Standard Model processes such as 
the production and decay of the Z boson and lower mass resonances such as 
the ^/ip. An essential step for the discovery of new physics will be to prove 
and improve the understanding of the Standard Model at the LHC. 

Once the detectors and the backgrounds are understood, the search for 
new physics, the Higgs boson, supersymmetry and others (extra dimensions) 
will commence. If supersymmetry is at low mass as hinted by electroweak 
precision data, large signals are expected and many observables will be mea- 
sured. 

Disentangling the fundamental parameters from the measurements neces- 
sitates powerful tools such as SFitter and Fittino as observables depend on 
different combinations of them. Additionally, the parameter space has large 
dimensions, thus efficient ways of sampling are needed. 

The project SFitter has developed a framework for such an endeavor. 
Tools such as weighted Markov chains have been developed to sample the 
parameter space. Interfaces to the most precise theoretical predictions of 
observables have been built. A rigorous treatment of statistical, systematic 
and theoretical errors, correlated or not, is provided. 

Using SFitter, two projects have been described here. The Higgs sec- 
tor was studied in the Standard Model as well as in supersymmetry, using 



67 



only the expected measurements for a low mass Higgs. The supersymmetric 
parameters have been studied at the electroweak scale in the MSSM. The 
parameters have also been studied at the GUT scale in a top-down approach 
(mSUGRA). 

In addition to these studies, SFitter has also analyzed the impact of 
{g — 2)^ on the precision of the determination of supersymmetric parameters 
in a Les Houches project |196j . In |197] the LHC potential for the parameter 
determination of Split-Supersymmetry was studied. New SFitter projects are 
already on the way. From one, excerpts were already shown in this paper: 
a study of the extrapolation of the supersymmetric parameters from the 
electroweak scale to the high scale. Another study will be the implication of 
the supersymmetric parameter determination at the LHC on the prediction 
of the relic density. Here we also plan to include direct /indirect detection 
results. 
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Figure 5.1: Extrapolation of the 
gaugino mass parameters from the 
electroweak scale to the GUT scale 
using LHC and ILC measurements. 



Figure 5.2: Extrapolation of the 
scalar mass parameters from the 
electroweak scale to the GUT scale 
using LHC and ILC measurements. 



The tools developed by SFitter can also be used for other models where it 
is difficult to disentangle for the measurements the fundamental parameters. 
Nature be supersymmetric or not, with the LHC/ATLAS now and an ILC in 
the future, particle physics will continue to provide exciting measurements 
and insights into nature. Maybe SFitter will take us, supersymmetry existing, 
with the LHC and ILC to the GUT scale, close to the Planck scale, as shown 
in Figures 15.11 for the gaugino mass parameters and in Figure 15.21 for the 
scalar mass parameters of the first and second generation. 
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